‘or lts. oe 


hon on Vb ivind ae 
Va Hi HS Ba Ver ie teeat a rset qhaiet 
{ ,' al + ‘ ‘ i) wv ee , : pp te S te 
i 1 es why 2 i i A Ny ri POC RR Hebe bid tharos ce 
‘ i { Rati t VERE Rein aycg Cey Hb) Rey johabing ct 
i Wily dot a ' HPA C a Nobaeoqtermends 4 Ayicnk Prades cyte ase MARIANAS bale me 
i ita ; fatite sae , 2Laheyo) Seghina cewee thy 
‘ 4 \ 4 s ¥ ) Un WS ¥ ‘ bee i te a Ces pe) rt By ee 
; ‘ hg CPR OV AAS OS, i Y ‘ 
‘ f i ' 4 4 tay i uy i 
' i "4 satay i] " i 4 
4 } 7 it 1 Ky ' ’ ere 4% iW 
1 ’ ‘ ' \ " rn ’ ‘ is 
3 ‘ \ « ‘ i rT ’ ' s “ 
; H ; : ; AS ast Ey iat Steins hin¥ ul yes 
AS re Pyeng i aay ee FACASTAN A SoA OwanuaR eens 
; icy : : HACE AL core Saeeey testa gets 
‘ v4 rae) ayy RGR < Funes oe 
: OA yr ineny Hi Ne $i SET yeNrhos by 
eet : ‘ al st ; ys WHEATON TS heli sR WaAA SAE eLS 
f ip ‘ aS i ; We bah BS Coal anes, 
i ’ ' aR Tay Ent i i 
a pets BTR ME Stee Oita saree 
: i { rt ‘ ’ rn i Pye 
; : TA ' y Citra kes tadive ge bai baala 
frien ; ; ial i mahes WAG Has g 


For Reference 


NOT TO BE TAKEN FROM THIS ROOM 


; 6 j reer Spine Met eb : 
: ' : ‘ 
4 1! i A a *Y A f, ‘ AWanrey ‘ 
i I : regs Pih.8 $ } LAT OT od at trie 
t ‘ , biased bixeesed i teil a4 c8N eAp Ads v0 
} I bs int sd hans Da pectieed ty igsned irae swe 
{ we vi i i eke art 
i ‘ eitai f ve . Aerieend WAL RE Ope: 
8 ! j wii “ op teed hee. 
: , : Ms diets itbrce 
1 oid ba dal , t {Oise Ades aoe ce ‘Reet, 
} fe { § a rane 
i i bobs tayshan dey MVE *8 0a Ld Geuodanrgonceg is 
; lb ‘ i Fi rE Meee} bei kidaen Dio tte tee ie 
‘ s ‘ Pere : J rags talgnd soot Etta ehehe Gt idlueh iteheipetc ly 
} ! 466 : M4 ips RRS igitiredens Sanehsegnantee tere ta" 
; y , if dite } { Nandaitete HUES LURE VET ae e rN a aio Gu re ehegeae ied s 
{04 ryihg Uae CHE h LAMA EL hs 48sn alas staiad eRthavsbebetiee fs 
i i ; iS t Bde 4 Maeva tus nas ce aig Vedtacdthat sta ete Sight e 
c jee’ Kde ae hah yd uae pfudce keg ans ye es MOR eas 4474+ eQsicoeinaeg Ierighoemion 
: i i ‘ ACSA Aan ate neg dt Py eh-heueev renin: ees fey desdewen hits ¥ets+G0) 
bh key is Lea td awed sh 443 Wolsstieeni seit fUtieooe es. 
; +}: Ry a ' br desketey aTytt hcaegst Gok Seae DP RE aD INA bere 24anarsh 4 
f pats Hiatal ate ated th deates MOREA te Foes PR ais oS tort en sieh 4 
2 tas iy i CEH OAD Pa eae BUND Bt £1) perry WIth sFrhslavegeesdee 
t ; Ph tc Ai Hess et Rubay ease wabaHtih rhaloust caus aoe cup eet 
mh ) i he SOR hadedate tyne nvagh Vine tatpiniad ha tenn tub eked 
P fie: ‘ eg Veit bdeee a PARR Teeth: 
i jst ‘ ‘ heat Ciel ed ecoeey Sica theks reaction eben cee eRd 
j A AES sY YSERA Te ay nyt rhe’ Pabatee Rete en ae 
: : \ ap mega ea, fig § Mehbioseratiria tral wlaveeng aie 
hed ide! ie don WUWeA phate heen g bees gi ntsCaeckchotel hey sheet 
hy ref i i ‘ bv yrk eye WORSEabe8 ¥-b¥ 954 
{ i { ' i heat PST 11 sheledatou rabies 
{ ! ‘ MiG ae 
deh cbse Hagid alad Tobe ; 
‘ Y : iw; FSS ee th ore ee wart, Wonnartetagh Dida be COPD PP ene ae am 


Gx wenis 
UNIOEASTTATIS 
MBERTAEDSIS 


THE UNIVERSITY OF ALBERTA 


RELEASE FORM 


NAME OF AUTHOR Bruno SAM YUE CHI 
TITLE OF THESIS MOMENT REDISTRIBUTION AND SECONDARY 
MOMENTS IN TWO-SPAN PRESTRESSED CONCRETE 
BEAMS 
DEGREE FOR WHICH THESIS WAS PRESENTED MASTER OF SCIENCE 
YEAR THIS DEGREE GRANTED FALL, 1984 
Permission is hereby granted to THE UNIVERSITY OF 
ALBERTA LIBRARY to reproduce single copies of this 
thesis and to lend or sell such copies for private, 
scholarly or scientific research purposes only. 
The author reserves other publication rights, and 
neither the thesis nor extensive extracts from it may 
be printed or otherwise reproduced without the author's 


written permission. 


THE UNIVERSITY OF ALBERTA 


MOMENT REDISTRIBUTION AND SECONDARY MOMENTS IN TWO-SPAN 
PRESTRESSED CONCRETE BEAMS 
by 
Bruno SAM YUE CHI 


AMtHESIS 
SUBMITTED TO THE FACULTY OF GRADUATE STUDIES AND RESEARCH 
IN PARTIAL FULFILMENT OF THE REQUIREMENTS FOR THE DEGREE 


OF MASTER OF SCIENCE 


DEPARTMENT OF CIVIL ENGINEERING 


EDMONTON, ALBERTA 


FALL, 1984 


< ae ey a 
opae te et a 
ary oP Ata, 
y Net oh SE ee Or 
1 en Wee 7 
7 : - 
ATREEIA IO YTESHAVIMU BHT 
; - on : 
“ATZ-OWT WI 2THEMOM saxanossy ama woemaiame tase 7 
eitacaeliadiaaictiad hpassateas | | | 
Ne eh a, 
THD BUY MA@ one oo At 
) Fas - Pes, | PAL ba 
z : 4 we r is” 7 
ft © 
i : i as: 
atom Hero gl ‘ae 
HORASZAA twa a igure statiasto #0 Ya sits tty a aareu a = 


ely 


‘waaouc ant TO area 10peE ‘get +0 TEMS eat ahctanl 
rs save toa ‘3G. ABTBAM 70 


2 ar OF rs ® as ¢q 
7 3 = van i ty ic _ i : vA i : af ks fs tics 
oti outset? ou aver <2 Sema ream 
: : ae, a re a <u 2 ‘gi a de av ee . > ‘oe a ed 
i Sf © tay | e = ip 7 ?. 


THE UNIVERSITY OF ALBERTA 


FACULTY OF GRADUATE STUDIES AND RESEARCH 


The undersigned certify that they have read, and 
recommend to the Faculty of Graduate Studies and Research, 
for acceptance, a thesis entitled MOMENT REDISTRIBUTION AND 
SECONDARY MOMENTS IN TWO-SPAN PRESTRESSED CONCRETE BEAMS 
Submitted by Bruno SAM YUE CHI in partial fulfilment of the 


requirements for the degree of MASTER OF SCIENCE. 


ABSTRACT 
The object of this study is to investigate the flexural 
behavior characteristics of continuous prestressed concrete 
beams and to determine the effect of inelastic behavior on 
the secondary moment. The analysis uses conventional 
moment-curvature relationships and compatibility of geometry 
in predicting the complete moment-load curve to failure for 
a given beam. The analysis was compared with experimental 
results. It was found that a variation of secondary moment 
is likely once cracking occurs. Finally, from this 
investigation a design proposal involving the magnitude of 
the secondary moment at the ultimate limit state is 


recommended. 
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1. INTRODUCTION 


1.1 General Remarks 

The subject of moment redistribution in continuous 
prestressed concrete beams, and its effect on the secondary 
moment produced by prestressing, has become a matter of much 
discussion. 

A variety of approaches has been used to demonstrate 
that the secondary moment should be considered when 
calculating the ultimate load capacity of the beam, and can 
be neglected if and only if full redistribution of moment is 
achieved. Complete redistribution of moment is not likely in 
most practical cases. There is no experimental evidence that 
the full inclusion of the secondary moment in such instances 


yields a safe design. 


1.2 Object and Scope 

The purpose of this research is to study the flexural 
behavior characteristics of continuous prestressed concrete 
beams from the post-cracking stage up to ultimate in order 
to establish the relationship existing between the inelastic 
behavior and the secondary moment. 

The procedure used is to set a theoretical model that 
is capable of tracing the post-cracking behavior of such 
beams, based on fundamental concepts outlined by previous 
investigators. The analysis described in Chapter 4 involves 


the use of moment-average curvature relationships and 


ai yietil somet tesmont % OFF MEET IEAOS: esolqmod .bsvetdos 
eds earabive Lednamiisaxs on ai evecdT 89289 Lesitoerg Jeom 


seonsreni, dove ai inemos a aaa. Pre % setevioni stud odd ic 


oo 
Ls, 
a) 


avaunizno>. ai clitiaieala aremon 19 sala eat ae | | 
yisenosse edz no soetis eps bone aneod stetsnea bee: = | 
daum 20 tedYenm « satoned eur (eoiganizes3s Ve besuboxq : 


snole 


eseticnemeab of Hseu need aay peso legge, i090 yipsasy A 


.7 a. _ 


natiw batshleqsa sd Bivors tasmom qisbno2e8 ad? 36 _ 
en bine aod aie jo qhiosgsd beot. siewistu ods rol salvation . 
ai tasmom %6 soisudiage thes tit ii. vite bre 2! bes3elpen | 


7 


7 
a) 
a 


= ne 


~~ 


a 


ts 


neieeb s2ne 5 


~ 


ages? bas soetee a ie 


Nak 
7 


} ) Ws - 
feruyst? Sis ybere 3 ee dpxseas1 efde $0. eaegsg _ es _<—& 
an i 


ststonop besastdesiq: auounisne> in edi se inganesasts 10 


WwHIo. ni a aa ar ghitieeis-s 


principles of geometry to determine the distribution of 
moment at any location in the beams at various stages of 
loading. A check on the accuracy of the theory is performed 
in Chapter 5 by comparison with available test data. Design 


recommendations are suggested as a result of the above 


analysis. 
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2. DEFINITIONS AND PROBLEM STATEMENT 


2.1 Statically indeterminate ao netecctipT 
Although most prestressed concrete construction at the 
present time consists of statically determinate beams and 
girders, there are important advantages associated with 
indeterminate structures of prestressed concrete: 
~ Design moments are smaller for given spans and 
loads than for determinate structures; 
. Stiffness is increased and deflection is 
reduced; 
2 By continuing post-tensioning tendons over 
several spans, fewer anchorages are required; 
Joint rigidity available in continuous frames is 
an important mechanism to resist horizontal 
loads such as are induced by wind, or seismic 
forces; 
= Many ingenious arrangements have been developed 
to avoid the high frictional losses of 
prestress. 
As a result of these advantages, the applications of 
continuous prestressed construction are expanding, and this 
trend may be expected to continue (Ref. 1, 16, 17). Two-way, 
continuous flat plate slabs are widely used, and have proven 
both functional and economical. For medium and long span 
bridges, the economic and esthetic advantages of continuity 


are dominant considerations. 
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2.2 Secondary Moments 

When an eccentric Rrestressing force iS applied to a 
statically determinate beam as shown in Fig. 2.1, bending 
moments P.@ are induced; e€ is the eccentricity of the 
resultant tendon force P with respect to the centroid of the 
cross-section. The beam will deflect when prestressed, 
usually cambering upward, but no external reactions are 
produced by the prestressing force. 

For a statically indeterminate beam as shown in 
Fig. 2.2, the action is more complex. The moment just 
described, which will be referred to as the primary moment, 
induces a deflection as before, but the beam is restrained 
by the redundant system of Supports. Reactions are produced 
at those supports, giving rise to secondary moments in the 
beam. In this case, the total moments produced at any 
section by prestressing is the sum of the primary and the 
secondary moments. 

The magnitude of the secondary moments in any given 
case depends on the particular tendon profile selected. For 
special cases such as the concordant tendon' case, the 
secondary moments may be zero. They are usually comparable 
to the primary moments and in many cases may be larger, even 
though they are called secondary. 


' When the tendon profile selected produces no reactions due 
to prestressing, no secondary moments are developed. The 
thrust line produced by prestressing coincides with the 
steel centroid line, as would be the case for a single span, 
statically determinate beam. Such a tendon is called a 
concordant tendon. 
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Concrete centroid 


Steel centroid 
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Deflection due to prestressing 


Figure 2.1 Statically determinate beam 
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Concrete centroid 


Beam profile 


Steel centroid 


Se Deflection if centre 
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Reactions at supports 
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oe ee Bee Serene Actual deflection 
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Figure 2.2 Statically indeterminate beam 
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2.3 Treatment of Secondary Moments due to prestressing in 
the ACI Building Code 

The proper treatment of secondary moments in the 
ultimate load analysis has been the subject of much debate 
(Ret te18F u21, 6°22) 

In the 1971 edition of the ACI Building Code, it was 
stated that the effects of moments due to prestressing, 
including secondary moments, shall be neglected when 
calculating the moments corresponding to factored loads. It 
was also stated that the behavior shall be determined by an 
elastic analysis, with only a modest amount of 
redistribution of moments due to plastic behavior permitted. 
The accompanying ACI Code Commentary stated that the 
secondary moments produced by the prestressed force in a 
non-concordant tendon disappear at the capacity at which, 
because of plastic hinge formation, the structure becomes 
statically determinate. 

The 1977 Code, however, requires the consideration of 
secondary moments, using a load factor of 1.0, up to and 


including the ultimate load. 


2.4 Problem statement 

The contradiction arising from the Code consideration 
of secondary moments in ultimate state calculations may be 
misleading to the designer. Secondary moments need 
definitely be included in the elastic analysis. However, the 


load carrying capacity of a continuous member is not 


= . | iit 


_ 


eft al f#inemom cusbaoata! 39 ssensenss ogee | 
sjedeb dou bo t2etdue eile aged ger cieylans beet * ul 
Ss 8 at ‘The 

aoe 3% yShed Ere Bliee ISA one Io noizidbse rer on3 ee t 
eT A, 


enigusyeesig 02 aup eanainee 5, ag7s7 6 oda tan? | saab 


ngdw bese sipen set fedg: seaagmem enabs0¢8 entt nike | on 
ote 


a l 

] 
al : 
7 


31 -absol bs10td= sy oJ pithnanpettes ainsdon oft prizel 
6 yd _benimretsb ad. [Lerte IS£V5 tiga “ae? gary beset oats | nas hs 
2s sesteins fastion, se yins atie j2seylens si semds 

Ser1 ime voiwadad digesta ad, J arnentos to oottudixs 
eta 3sc batave Tyst eommod sho? 3A ent yRequemnney 

s nt. 60761 poaeersassg ont we oe ainemom Y18 
“ised 75 (2p0koRD ed? a TAIGES TP aba gt jnedt0: 
aemo ted sud 130 22 on oped nine’ apoid oigeeig io ; 
: te Pe ae, stentwuszeb yileg 3838 

46 ng? sezebienen bet soskypey favausd 8500 eet —_ 
bre 6s au went io ea382 patna 6 eniau Sanam is . 


7 ? 
- 
= 


affected by the secondary moment if complete moment 
redistribution can take place at ultimate. If plastic hinges 
do not fully develop, then ultimate load capacity will lie 
between the load resulting from an elastic analysis and the 
load at full redistribution (Lin and Thornton, ref. 18). 

As required by the present Code, secondary moments are 
to be considered up to and including the ultimate load. 
Since certain circumstances exist when some redistribution 
of moments does occur and is, in fact, allowed in design, do 
the secondary moments vary under these conditions? The 
question arises due to the fact that, depending on the 
System of supports and the tendon profile, the full 
inclusion of secondary moments at ultimate may lead to: 

> uneconomical design in sections where the 
secondary moment is unfavorable to the nominal 
Capacity; 

= unsafe design because the contribution of 
secondary moments to the beam capacity may have 


been improperly assumed to exist. 


It is worth noting that, from a design standpoint, the 
ultimate stage does not necessarily correspond to the 
formation of a mechanism, but rather to the loading at which 


ultimate capacity has been reached at a critical section. 
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3. LITERATURE REVIEW 


3.1 Nonlinear analysis of continuous beams 

The 1971 ACI Building Code provisions concerning moment 
redistribution in continuous prestressed beams are 
definitely inconsistent (Lin and Thornton, 1972, ref. 18). 
According to the Code, full moment redistribution at 
ultimate is not permitted, while secondary moments must be 
neglected at the same time. By means of examples, it is 
demonstrated that neglecting secondary moments may yield a 
non-conservative result. A method for determining the 
ultimate load capacity of a continuous beam is proposed. The 
method takes into account the secondary moments without 
calculating for them and the final moment configuration is 
an intermediate stage between the elastic case and the full 
redistribution case. The method is conservative due to lack 
of analytical and experimental research concerning the 
plastic behavior of prestressed concrete beams with 
non-concordant cables. 

It is believed that an exact solution can only be 
obtained when the moment-curvature relation for the entire 
beam is analyzed beyond the elastic range and up to failure. 

A series of tests of seven simple-span beams and three 
beams continuous over two spans of 28 ft each has been 
conducted (Mattock, 1971, ref. 21). Although the primary 
variable of the study was the effect of bond on the behavior 


of post-tensioned concrete beams, a considerable amount of 
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redistribution of support moment has been observed at 
ultimate. A large portion of redistribution (up to 85 
percent) has been attributed to the action of the 
non-concordant tendon, because the test beams had a net 
reinforcement index (w + wp) - w') that did not allow any 
adjustment of support design moments according to any 
edition of the ACI Code. Therefore the secondary moments 
have been assumed to have a direct effect on the amount of 
redistribution available. It was concluded that for a 
downwards transformed tendon profile, "redistribution of 
design support ultimate moments by an amount equal to the 
positive secondary prestress moment should be allowed in 
design, without a special limitation on the amount of 
reinforcement". This reduction in support moment does not 
require any inelastic deformation at the support section. 

These findings led to subsequent changes to the 1977 
ACI Code, which required the inclusion of secondary moments, 
using a load factor of 1.0, up to and including the ultimate 
state. 

Many theoretical approaches were developed to enable 
the distribution of moments at ultimate to be related to the 
physical properties of the beams and the pattern of loading. 
These ranged from Guyon's general analysis which takes into 
account the actual distribution of curvature along the 
length of the beam (1960, ref. 14), to Baker's simplified 
approach in which the inelastic deformation is considered 


concentrated at the critical sections, i.e. the concept of a 
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"plastic hinge" theory with the "hinges" having limited 
rotational capacity. 

It is shown that yield of reinforcement can provide 
advantageous moment redistribution. However, the small 
amount of steel yield available in prestressed concrete 
beams may reduce the moment redistribution possible as 
compared with ordinary reinforced concrete (Baker, 1949, 
nefia 2); 

In general, the theories require a knowledge of the 
moment-curvature relationships for the beam sections. The 
theory developed by Priestley et al. (1971, ref. 24) takes 
into consideration the variation of curvature between cracks 
caused by concrete tension. It also showed close agreement 
with experimental data. The relationships between moment and 
average curvature have been used to determine the 
moment-load curves for continuous beams up to the onset of 


concrete crushing (Priestley and Park, 1972, ref. 25). 


3.2 Rotational capacity of hinging regions in reinforced 
concrete beams 

Instances can occur in which the strain capacity of a 
reinforced concrete hinging section is exhausted before full 
redistribution of bending moments is achieved in the 
structure as a whole. It is therefore necessary to consider 
the deformation of the hinging regions in any theory of 
limit design for structural concrete, and more specifically 


to limit their rotation to known safe values (Mattock, 1964, 
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The main factors relating to rotations are: moment 
gradient, concrete strength, reinforcement yield stress, 
beam effective depth, amount of tension reinforcement and 
confinement of the concrete in compression (Corley, 1966, 
ref. 12; Roy, 1964, ref. 26). 

It is demonstrated that similar approaches could be 
used for both reinforced concrete and partially prestressed 
concrete in evaluating moments and curvatures (Bishara and 
Brar, 1974, ref. 3). 

Computer-Simulated flexural tests carried out to 
identify all major variables that affect the behavior of 
partially prestressed concrete sections, confirmed earlier 
findings concerning the ductility of ordinary reinforced 
concrete sections. For prestressed sections, ductility 
showed considerable sensitivity to the effective 
prestressing. However, variables such as cross-section shape 
and high-grade steel stress-strain relationship have a 
relatively minor effect on the inelastic behavior (Cohn, 


ioG2zee cote. 11). 


3.3 Experimental programmes 

Only few experimental data were available in the 
literature. Major extensive test programmes on prestressed 
concrete beams were conducted by Warwaruk on simply 
Supported beams (1962, ref. 28) and by Hawkins on two-span 


continuous beams (1964, ref. 15). In general, flexural 
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cracking and bond were major factors in the behavior of test 
beams. It was observed that moment redistribution in 
continuous beams is initiated by flexural cracking and 
relative reduction of stiffness over the interior support. 
Pronounced redistribution occurs only after the moment over 
the interior support reaches the nearly flat portion of its 
moment-curvature relationship. The development of inclined 
tension cracks reduced both the load carrying capacity and 
the ductility of test beams failing in shear. The basic 
mechanisms of failure in shear or in flexure were similar to 


those in simply-supported beams. 
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4. DESCRIPTION OF THEORY 


4.1 General remarks 

In order to obtain a better understanding of the 
behavior characteristics of continuous prestressed concrete 
beams subjected to bending, a theoretical model has been 
developed to study the phenomena occurring in the beams when 
loaded up to failure. 

The theory requires a knowledge of the moment-curvature 
relationships for the sections to determine the distribution 
of moments throughout the structure at a specified load. The 
moment-curvature relationship formulated by Priestley has 
been adopted. The theory takes into consideration the 
variation of curvature between cracks caused by concrete 
tension and makes possible the prediction of both the 
curvature at a crack and the average curvature along the 
length of the member. 

A computer program using Fortran statements and based 
on the above theory has been written. Further details of the 
computer program are enclosed in the Appendix. Comparison 
with experimental data is presented in the following 


Chapter. 


14 


efi do paloqetessbay ery & niasdo, co: 
eis7ono2 haagetIea7G sxeumigao? 30 aah atanont 


neod ent Iebom feod seicet? ‘ saniBoon, ce 
siener of | 


) .souLiat e2 qu Fk 
to ab 35Wo be * astiypss sweet set 


netqudisath sas snimiests> on 2 Se do ad? 202 aq teh ane 


ea7 .bsol. bei?rooge 6 ts magus a: ate gues 


3 aiagaotseler gu 788 


ngdiw emoed a3 of geitaueces Bnaqisaee etd, 


Me) 


stu 2 eY wo DS ASCs eng 


en yoltéeir® yas Beatie Lu 
tsesebizrios ce eae groans — +b 


sos oy ties oor, alias ean ea 
er3 aged 16 CLI 9039) odd TENA ETA: sa ad 


Tut ope: rove aii has. set 6 3% 


bensd. bus 234% aipte Hero priek Fars en ris 

A or 7 " 3 am =< ; - 7, 

enc toed teat ter) 103 -ngiatam palin “eat ore F = 
- : Ts ais ¥ 

nde tragmay:” wahnan ‘beeot 4s ie 


“ wolie? siz ft beineaoig es sent 
a 


aii pacis ¢ 


4.2 Analytical model 


4.2.1 Assumptions 


1) 
2) 
3) 


4) 


5) 


6) 


7) 


8 ) 


Quasi-static loading; 

Bonded beams; 

Negligible shear effects; 

Plane sections remain plane (linear strain 
distribution); 

Any known material stress-strain relationships; 
Effective (after losses) prestressing; 
Linear-elastic behavior up to decompression of 
concrete; 

Partially prestressed beams in which limited 
cracking is permitted by the designer at service 


loads. 


4.3 Moment-curvature relationships 


4.3.1 Prestressing steel strain and stresses 


In prestressed concrete beams, the prestressing steel 


Strain is not zero even though no external load has been 


applied. It is therefore necessary to include the initial 


strain in the analysis. 


Strains in the concrete and steel at loading stages of 


imterest are shown in/Fig. 4.1. Strain distribution (1) of 


Fics. teresults from the application of prestress force P. 


acting alone. At this stage the stress in the steel and the 
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Figure 4.1 Strains in concrete and steel 
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associated strain are, respectively, 


ot 

f, = — (aie 
A> 
: ae 

€; = = (4%.2)) 
Fa 


The steel strain is shown with respect to its own separate 
Origin. Stage (2) corresponds to the decompression of 
concrete at the level of the steel centroid. Assuming that 
bond remains intact between the concrete and steel, the 
increase in steel strain produced as loads pass from stage 
(1) to stage (2) is the same as the decrease in concrete 
Strain at that level of the beam. It is given by the 


expression: 


€2 = re ( + =) (4.3) 
NewS 


in which g is the radius of gyration of the cross-section. 


When the load is increased further to the stage (3), the 
neutral axis is at a distance Cc below the top of the beam. 


The increment of strain is: 


em Die. 24 (4.4) 
e 


The total strain is the sum of the three components 


a he haath Sy: +N Go TEs (4.5) 
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and the corresponding steel stress is given by the 


stress-strain relationship of the particular steel grade. 


4.3.2 Cracking point 

One loading stage that is of interest in a section 
analysis is the one at which cracking moment is reached. A 
Significant change in slope can be observed at that 
particular loading stage in a typical moment-curvature curve 
as shown in Fig. 4.2. The model will therefore consider the 
cracking moment as a benchmark between elastic and inelastic 


behavior. 


4.3.3 Conditions just before cracking 

Assuming that the concrete has a tensile strength, the 
conditions just before cracking can be defined as the stage 
at which the extreme concrete fibre in tension reaches the 
flexural tensile strength f,', so that cracking would 
certainly occur should the load be increased by a slight 
amount. 

For a rectangular section as shown in Fig. 4.3, 

Concrete compressive force 
Cc 
C=b f.edy (4.6) 


Cc 


Concrete tensile force 


fan (ieee Gob (4.7) 
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Figure 4.2 Moment-curvature relationships at a section ina 


constant-moment zone (Ref. 24) 
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Figure 4.3 Conditions just before cracking 
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Steel tensile force 


The concrete stress-strain relationship is given by the 


expression of f in terms of e. The equilibrium equation is: 
Gee Tectidp (4.9) 


After each component of the equation is substituted 
with an expression involving only the neutral axis Cc, and 
the steel strain ep, the equilibrium equation is solved 
Simultaneously along with the steel strain equation 4.5. 


The curvature is derived from the strain, distribution 


{ EL, +e. 
({) . aires CATO) 


and the corresponding moment 


Cc 


M = b/f ydy +b £,' S(n-e)? teApfp, (deo) = (4.171) 


O 


4.3.4 Conditions after cracking 
When the stress in concrete extreme fibre in tension 
has exceeded the tensile strength, cracking occurs in the 


section. As a result, the component 7, no longer exists. 


qidaroitelet atexse~ -anetis ot 
sit.» to @m7s? ai ie soit 


prt wily. 8 
e453 yd neve 2t 


ei aotteupe mu isdii beps 


a - SP 8 se 


(eb) gt E 


to Jieguqtos daae sn 


| . - Ad; 
Ketotreedue @i- MOLIeyps eee: 


ET enews fayn? |G? senvaat Ha 


yw 


2 ea: notisups aprsastinle ott . qo bande le 


: ; : . se ot ok tsi & 
cuns- t¥Yetze Inege Sad aeew ae) aoe bps 


. ) . j pavdasd ‘ef 2% pw iT 
‘" ade ek pheene sag me1] bevegeqd 2: 2 rt ey o] 
ta’ a het od a a he, i AD 
, | £5 
‘ A ; 
or. 7’ f 
pi & ing % > Saal ; 
“40 > — ee = ny 
4) 
{ : z a 
: 1 4 
2 7 
4 : i 7 


1 Acts teenie onteriogge1%02- : 


) : 

5 
6 . 
a . pes , ee 


J ‘ a a rina _ % ; oe ae. 


7. 


oe. : m, an a A. yok } ‘a 
4s 1.2) a=) yg 4h ¥ Mae ne fe * VRE z as 


‘. 


on . 1a : ' 5 Sead ain resis sensi 
eal notanas perres snosians wos 
si 


aici ae 


ne ciel nzerse = ot " 


_ gy 


22 


The location c of the neutral axis can be determined 
from the general cracked section analysis developed by 
K. Shushkewich (Ref. 27). The equation of neutral axis is: 


zbNe* + SbMe? + (BN+aM)c - (yN+BM) = 0 (4.12) 


in which the different coefficients can be calculated from 


the general transformed section shown in Fig. 4.4: 


a = (b-b,)h, + npAp (4.13) 
B ee ie be NpApd (4.14) 
Y = L(b-b,)hy + npApd? (4.15) 


where b,, = b for a rectangular section. 


N = F = ApEp (e,1 + €2) (4.16) 
M = Mé®-eFd (4.17) 
The stresses are: 
Concrete 2 a eee (4.18) 
6 
Steel Ep gies ee 2 (4.19) 
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Figure 4.4 General cracked section analysis 
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4.3.5 Average curvature 


4.3.5.1 Bond stress 
It 1S generally agreed that the concrete-steel bond 
stress iS a maximum very close to the crack and 
decreases in some fashion further away from the crack. 
For this analysis, it will be assumed that, at 
first cracking, the bond stress decreases linearly from 
a maximum U, at the crack to zero at a distance away 


from the crack as shown in Fig. 4.5. 


4.3.5.2 Bond length 

Immediately after the first crack forms, when 
M=M._., in a region of constant moment, a stress 
condition that varies between two limits exists. 

Let A be the section at a crack, and B be a section 
some distance away from the crack where the stresses 
have not been affected by the formation of the crack. 
The stresses and the stress resultants at section B are 
as described by equations 4.6 to 4.11. At section A the 
stresses are found from equations 4.12 to 4.19 with 
M = M., - 
At section A, all the tensile stress is carried by 
the steel; between sections A and B, tension is 
transferred from the steel to the concrete by bond. If 
f.q and f,p are respectively the steel stresses at 
sections A and B, then the minimum distance #p from the 


crack over which sufficient tension can be transferred 
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Concrete 


Cracks 


Prestressing steel 


Figure 4.5 Bond stress distribution between adjacent cracks 
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from the steel to the concrete by bond, to cause the 


modulus of rupture f,' to be just reached at B is: 


“ Ap (fa - fi) 


Lp % 
ee 


(4.20) 
in which Ugy is the average bond stress between the 
steel and the concrete, XO is the total surface 
available for bonding of the prestressing steel per unit 
length, Ap is the total area of prestressing steel. The 


maximum bond stress is: 


Um = 2 Uay = 2p (Fa fa) (45,.2'1 ) 


g, X90 


4.3.5.3 Crack spacing 

A new crack cannot develop between two existing 
cracks which formed when M = M... if the spacing between 
these cracks are smaller than 24p. Sufficient length is 
required each side of the potential crack position to 
build up enough concrete tension and induce a new crack. 
It is evident that, with initially random cracking, the 
individual crack spacing soon after the formation of the 
first crack will vary between two limits, lp and 24p. 
The average crack spacing will be approximately equal to 
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4.3.5.4 Stress distribution at a distance ¢ from a crack 
ist, 

At a distance £ from the crack, the steel stress 
will be reduced by bond and the concrete tension will 
build up. The reduction of steel tension force over the 


length & from the crack is: 


Js 
AF = [. ZO dé (4.22) 


O 
Therefore, if fs is the tensile steel stress at the 
crack for the particular moment M acting, the tensile 


stress in the steel at distance #& from the crack is: 


fe te ee (4.23) 
A, 
l 
oe ean kites 7 (4.24) 
Ap 
O 


From Fig. 4.5, the bond stress distribution can be 


expressed as: 


Up ==, ( = 7 | (4,25) 


Substituting u from eq. 4.25 into eq. 4.24 and 


integrating: 
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Substituting Um from eq. 4.21 into eq. 4.26 yields: 


2 
Pemeen ator oh) eles CAni2 7) 
Peete 
Note that f, from eq. 4.27 is independent of the 
maximum bond stress Um and thus the magnitude of the 


maximum bond stress does not affect the moment-curvature 


Curves. 


4.3.5.5 Average curvature 

The assumption is made that, at a section some 
distance from a crack, the steel strain is still 
linearly related to the concrete compressive strain. 

The known steel stress f, allows the determination 
of the conditions in a section at a distance & < py from 
a crack. 

The action of the decompression force, along with 
the external moment, can be represented by a resultant 
force R applied with eccentricity Crop above the top of 
the section as shown in Fig. 4.6 (Nilson, ref. 23, 

p. 97). The portion of beam can then be analyzed as an 
ordinary reinforced concrete member subjected to an 


eccentric compression force. 


Equilibrium equation: 


R. a6. Co-2T eT. (4.28) 


The concrete tensile force 7. can be expressed as: 
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tn =. CB= The-eR (4.29) 
in which 
Coe ae b +(-2o=he) y (p-b.) (4.30) 
Se 7 c3 |Owe C £ w ° 
T, = Apf, (4.31) 
R = ApEp (e1 + €2) (4.32) 


The neutral axis a distance c from the top surface, for 
the equivalent homogeneous transformed section, can be 
found from the equilibrium condition that the moment of 
all internal forces about the line of action of R must 


be zero: 


(d+e top) Apt, + To (ese +c +t hs 


3) 
= £63 (5 +e (b,c) 


hy es (ob, (B+ oe) 


Ang RE (orn, (SE oe (45233) 


+ 


+ 


with b, = b for a rectangular section. 
Solving Equation 4.33 with f, obtained from 
Equation 4.27, the position of the neutral axis can be 


determined, and hence the concrete stress at the top 
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surface using Equation 4.18. The curvature at this 


section is given by: 
i\ E& +E 
—— : (4.34) 


€. and e, are the strain corresponding respectively to 
the concrete stress at the top surface, and ff,. 
Integration of the curvature over half the distance 


between the cracks gives the average curvature: 


075 A, 
Te Soy Ai zs ; 


4.4 Deformation compatibility requirements 

In the analysis of continuous beams, the distribution 
of moments must meet the deformation compatibility 
requirements. 

The model beam that has been adopted considers the 
availability of test data. This allows a check on the 
accuracy of the theory. A set of experiments being referred 
to is from Hawkins' work (ref. 15), which consists of a 
series of tests on two-sSpan symmetrically loaded continuous 
beams. The investigation will thus be limited to a two-span 


symmetrically loaded continuous beam. 


4.4.1 Elastic conditions 
In elastic analyses of indeterminate beams, it is 


necessary that the slope at any interior support be 
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continuous. Consider a continuous beam as shown in Fig. 4.7. 
Neglecting any differential settlement of any support, the 


rotation @,;, occurring at j in span ij can be expressed as: 


—j 


O55 = y (x-x; )dx (4.36) 


x 


{ 
Lij 

on 
where the curvatures (1/r), are evaluated at each beam 
section from the corresponding values of bending moment, 
uSing the moment-curvature relationships defined in the 
uncracked section analysis. Similarly, the rotation 6;, 
occurring at j in span jk is: 


Xk 


6. = 
jk rc 


zi a (x, -x) dx (4.37) 
ie 
i 


Taking all rotations positive counterclockwise, the 
compatibility condition for the interior support Jj in 


Bia.e4.8° is: 
0 5 i - 63, = 0 (4.38) 


In the particular case of a two-span symmetrical beam loaded 
symmetrically, the compatibility equation reduces to: 


ji 


SG 
= Shixa=X pax =U (4.39) 
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Figure 4.7 Rotation of beam within a span 
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Figure 4.8 Rotation of beam at internal support 
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4.4.2 Non-linear conditions 

Another way of formulating the compatibility equation 
is to ensure that the deformations occuring within a span 
(between two supports) are geometrically compatible. These 
equations are described by Guyon. Consider a symmetrical 
two-span continuous beam with the load applied symmetrically 
as shown in Fig. 4.9. The beam can be assumed to take ona 
polygonal shape at an advanced stage of loading. The 
rotation at support B remains zero due to symmetry, but a 
Slope @' occurs at some small distance m from B. The final 
bending moment diagram, after the iterative process 
described in Section 4.5, is such that the deformations 
resulting from it are geometrically compatible. A set of 
bending moments can be considered correct when it satisfies 


the equation 


o(rA£) = 6' (4-m) (4.40) 


The quantity (4-m) will be approximated to 4 because of the 
difficulty to evaluate exactly the distance m, which varies 
with the width of the support, among other factors. Equation 


(4.40) is then modified to 


¢ A = 6' (4.41) 
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Figure 4.9 Compatibility of geometry 
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This compatibility equation provides a simple means of 
establishing an acceptable distribution of moments at a load 
higher than the load at first cracking. An error margin of 

+ 0.05 radians will be adopted as a result of the 


approximation described above. 


4.5 Procedure of analysis 
The successive steps of analysis are: 

1. Compute the cracking moments and the ultimate moments at 
all the critical sections, based on the section 
properties; 

2. From an elastic analysis, determine at which critical 
section cracking occurs first, and check the 
compatibility at that particular loading stage using 
equation 4.39, that is: 

a. From the resulting bending moment diagram, the 
corresponding curvature at the centre of each 
segment is found by referring to the 
moment-curvature relationships defined in 
Section 4.3.3; 

b. These curvature values are then used to determine 
the left-hand side of equation 4.39, which will 
equal zero if the assumed moment configuration is 
correct. Usually this will not be the case and 
adjustments are then made to the bending moment 
diagram by reducing either the maximum hogging 


moment value or the maximum sagging moment value 
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until the left-hand side of equation 4.39 does equal 
‘zero aS required; 

3. The ultimate limit state is established as the load 
corresponding to the final distribution of moments that 
Satisfies equation (4.41). The first iteration is 
performed using the ultimate moments calculated in 
Step 1, one of these moments being reduced in subsequent 
iterations. In addition, the moment-average curvature 
relationships defined in Section 4.3.5 are used in 
cracked regions; 

4, Stages between cracking and ultimate can be identified 
using intermediate values of moments. By modifying one 
of the moments and holding the other constant, 
compatibility can be satisfied, and the resulting 


distribution of bending moments can be calculated. 


One can then obtain a series of points on the 


moment-load relationships of the critical sections. 
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5. RESULTS OF ANALYSIS 


5.1 General remarks 


The moment-load curve has been established for a series 
of beams, based on the theory described in the previous 
chapter. 

In order to evaluate the accuracy of the analysis, some 
beam data have been taken from experimental work on 
continuous prestressed concrete beams, so that comparison of 


the results can be made. 


5.2 Experimental results 

Most of the experimental results available have been 
done on simply supported beams, due to the early research 
interests which focused mainly on the flexural strength and 
deformation characteristics, or the effects of bond. 

An extensive literature search provided only one 
detailed testing programme (Hawkins, 1964, ref. 15) on 
continuous prestressed concrete beams. Tests were carried 
out on 22 two-span continuous beams loaded at the midspans. 
Because the purpose of the investigation was to study the 
action of both bending and shear, with the emphasis on the 
effects of shear, the test beams were designed with varying 
amounts of shear reinforcement along with flexural 
reinforcement. In classifying the modes of failure for the 
test beams, the criterion used was the crack pattern 


observed. The six beams which failed in flexure were 
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selected for the present comparison. Further details on 
these beams are listed in Tables 5.1 through 5.3 and shown 
ry Gees Fe ge 

Tables 5.4 and 5.5 give a summary of results from both 
analysis and experiment. The only known values from the 


experiment are the cracking moment and the ultimate moment. 


5.3 Behavior of the test beams 

The position of the concentrated loads implies that the 
positive moment was five-sixths of the negative moment as 
determined by an elastic analysis. 

The first flexural crack was observed over the interior 
Support where the elastic bending moment was the largest. 
The appearance of this crack was accompanied by an 
adjustment in the relative magnitudes of the exterior and 
interior reactions. Beyond this stage, the moments at 
midspan and interior support deviated from the elastic 
distribution. Since the interior support section cracked 
first, its moment was gradually redistributed to the midspan 
section. When the midspan also started to crack, moments 
were redistributed back to the interior support section. As 
the load was increased the moment ratio remained essentially 
unchanged up to the theoretical ultimate load. The 
- moment-load curves shown in Figs. 5.2 through 5.7 include 
the effects of secondary moment as explained later in 


Section 5.4. 
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Figure 5.1 Reinforcement patterns 
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Figure 5.2 Moment-load curve for Beam 1 
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Figure 5.3 Moment-load curve for Beam 2 
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Positive Moment 


49 


hy GE “Fi 


1% ¥ 
: \ 
> Vas. . ‘ce 
NS VS 
a 
; s 
ae 
_ 
, 7 ’ 
oe z ; ; i 
. oF 
ch | 
4 
= eee 
4 ; 
AX 
) : a & re | 
‘tasaok vistegs = i ee 
SersmoNh SV AIS 22>, Net ek 


COE a Hes oar . % 7 oe 
ee ee eee Ree * 


a oS a os | oon 


ve a 


ei és ‘geass tot ad ae } 
ie eae ve 


4 ~ | fy ¥ 7 LS 
Gasnsadzeg a) pnbaey 
Oe . fone, rea 
" igiialatap) mee ‘oanataty 
7 Aesteleapesy sinks ve a8 hoses 
. te : = 


‘an iP, 


i esntuatio) 


Load,+40 Kips 


A A 
fe) 30 
8 7 
_ Wha 
< q 4! 
~ ve 20 7 
! ee 
\\ O° 
‘\ i 
x / 
SS 
\ 10 i 
Se vA 
Support Moment Be Ye Positive Moment 
=—300e =200°. —100) ~.\) 100 200 300 400 


: Kip-in 


Secondary Moment neglected 


—-—-— Linear distribution 
----— Inelastic behavior 
4 Ultimate Stage (experimental) 
© Ultimate Stage (calculated) 
- Intermediate value (calculated) 
+ First Cracking (calculated) 


Secondary Moment included 


Linear distribution 


Inelastic behavior 


Figure 5.5 Moment-load curve for Beam 4 
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Figure 5.6 Moment-load curve for Beam 5 
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Figure 5.7 Moment-load curve for Beam 6 
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Beams 1 and 2: No redistribution of moment occurred 
because the moment capacity provided at the critical section 
waS in accordance with the elastic distribution. However the 
tendon arrangement used in practice is likely to be 
different from the profile obtained by elastic moment 
analysis due to economical reasons. An example of practical 
design is shown in Fig. 5.8; 

Beams 3, 5 and 6: The moment capacity provided in these 
beams was equal at both midspan and support sections. Beyond 
cracking at the interior support section, the positive 
moment to negative moment ratio rapidly reached the value of 
unity. The correct distribution of moment was obtained using 
the compatibility equation 4.40. As the load was increased 
up to ultimate, the variation of moment indicated an 
inelastic behavior which was in reasonable agreement with 
the experimental values at the ultimate stage. Test values 
of moment and load at intermediate loads, however, were not 
available for comparison; 

Beam 4: A similar type of behavior was observed in beam 
4, in which the largest amount of redistribution occurred 
due to the tendon profile selected which differred 
substantially from the elastic case: the positive moment 
capacity to negative moment capacity ratio was 1.3, compared 
to the ratio of elastic moments of 0.83. 

In all six beams full redistribution of moments could 
be achieved as expected due to the low value of the 


reinforcement index ranging from 0.063 to 0.253, indicating 
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Figure 5.8 Basis for selecting tendon profile 
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underreinforced sections. However the effects of shear have 
to be considered because the development of inclined tension 
cracking may reduce both the load carrying capacity and the 
ductility of the test beams, as was observed in remaining 
beams of the test programme. The span-to-depth ratio of 9 
indicates the dominance of shear force which may result in 
deep beam action. Nevertheless, flexural failure is still 
possible. It appears that, with sufficient web 
reinforcement, full redistribution may be achieved even if 
the critical section or sections have very limited rotation 
capacity beyond yielding, as in beams 5 and 6. 

The higher values of moment and load observed at 
ultimate in the experiment may be attributed to the 
criterion adopted for failure: the ultimate moment 
calculated in the analysis used the concept of average 
stress of the concrete acting over the entire compressed 
concrete area above the neutral axis (Equation 36, ref. 28). 
The difference in tendon profile at the end portion of the 
beams between the analysis and the test might be one of the 
reasons for obtaining a higher load. However, this did not 
have much effect on the peak moment which always occurred at 
midspan and at the interior support due to the type of 
loading. 

Observations on test beams described in Ref. 15 show 
that the significant amount of redistribution can only be 
achieved through extensive cracking at the critical 


sections. Furthermore, the use of compatibility 
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equation 4.40 in the analysis has been made necessary in 
lieu of equation 4.39, due to the larger degree of 
deformation in the post-cracking range. 


“Sees 


5.4 Secondary moments 

As indicated in the moment-load diagrams, the amount of 
secondary moment existing in each of the beams was not large 
enough to be significant in the the over-all behavior and 
the load carrying capacity of the test beams. 

Secondary moments are calculated in Table 5.3 on the 
basis of the effective prestress force a . While the force 
in the tendon does increase significantly as loads on the 
structure are increased as a result of bending of the 
member, this does not represent a change in the prestressing 
force that produced the secondary moments. The force 
resulting from prestressing is unchanged, and the secondary 
moments are unchanged as the load increases up to the first 
cracking load. 

An intuitive approach can be established concerning the 
consideration of secondary moments: since they arise due to 
the restraint opposed by the redundant system of support, it 
is reasonable to expect a variation in the intensity of 
these moments as the structure exceeds the post-elastic 
range. There are several explanations for this: 

1) As cracking occurs at the critical sections, 
where peak moments are observed, the structure 


adjusts itself to the variation of stiffnesses. 
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This in turn produces a new arrangement of 
Support reactions. The portion of the reactions 
due to prestressing is likely to be different 
from the initial one because of the decrease in 
stiffness induced by extensive cracking in 
regions surrounding the critical sections. The 
secondary moments should therefore vary 
accordingly; 

As the loading proceeds further, one might 
expect the critical sections to form plastic 
hinges and reach the extreme stage of collapse 
mechanism. Secondary moments may then disappear. 
It is well established that they may be 
neglected at this particular stage (refs. 10, 
187022) 

In practical cases, full redistribution is not 
always possible because the combined action of 
bending and shear may cause one of the critical 
sections to fail by crushing of the concrete. A 
flexural failure, as usually defined, is not 
likely unless the shear-to-moment ratio is 


extremely small. 


| ie hd 


eois>oes eto 20 ne 


si pridoais svianesxe @e Babwon? = 


sit: .anoisoe? igsigiis 3a erie 


. o a 
“4 


tdpim ond ,tetragi 25a 99934 pribevt 

sttesin motos enerioes deoitegs act, rs 
adefion Io Bpsse 56 eT IKE ‘eke tne’ Sar euetain 

; c oe ea 

tnaeageei5S fists yam esnemor qrabnood@ -meinaizen a 

si Yom vedd Jede Ss det idade® iow et wo we 

) ; 7 3 

oaseen a 


Ot .#397) epsvé sEtvotsred Binh ieee 


= : 43s" i 
| | _ ~ & -+ * ; in ? xy 
r - 4 ' - 4+. mm os A. é ; ad i A ; 
@ le —— —S >t v ml - 
fon vt noriudintatbps 1lyg Skaeears eae 


$4 nsoaliso&’ Bacivnos sis szusgad sidt ae 


a 


[ssigris sds Jo PRO $28 4" gem eeee eas! 


- ‘ 


& ai o 2192 aris Grohl! 4 baie “octten 38 


ae a er 


ton 2k ibe! teb at id a L ve 


= - Cee 
»2@iL Ofgarounseton- 
ue? P ‘ee 


58 


For these reasons, the secondary moments may be assumed 
to decrease from the initial value, according to the 


following equation; 


3 


Pp = , 
cracking 
(Moc) t \ (Med i hee Pe ey ae 
post - crac full redist. cracking 


P is the load, usually the ultimate load, at which the 


secondary moment is evaluated. P, corresponds to the load 


racking 


at first cracking through an elastic analysis, and Pl vedic © 
the load calculated with ultimate moments at both critical 
sections. Also, the cubic variation is set to reflect the 
gradual reduction in member stiffness. Equation 5.1 is valid 
under the assumptions made when establishing the analysis. 
Thus the limitations are: 
1) (net reinforcement index) < 0.30 to ensure that 
ductility is available; 
2) bonded, partially prestressed concrete so that 
cracking is allowed. 
Equation 5.1 is graphically represented in Fig. 5.9, and is 
used in Figs. 5.4 through 5.7 to indicate the possible 
inelastic behavior when the secondary moment is taken into 
consideration. 
Caution has to be exercised when considering the 
yielding stage: although measurements of concrete strain at 
failure in beam tests indicate that values of e between 


0.003 and 0.004 are attained, a limiting strain of 0.003 for 
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Figure 5.9 Variation of secondary moment 
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the concrete will be assumed. Also, prestressing steels do 
not show a definite yield plateau. The spread between the 
nominal yield strength and the ultimate tensile strength is 
much smaller for prestressing steels than is the spread 
between the corresponding values for reinforcing steel. 
Therefore the small amount of steel yield available in 
prestressed concrete beams may reduce the moment 
redistribution possible as compared to conventional 
reinforced concrete. 

If the ultimate load does not correspond to full 
redistribution, i.e. if any of the critical sections has not 
attained its ultimate capacity, then a portion of the 
initial amount of secondary moment must be taken into 
account. This differs from the current Code provisions which 
require the inclusion of the total secondary moment. 

The behavior characteristics described in the present 
study have been demonstrated to be in reasonable agreement 
with the observed ones for a particular type of beam. The 
cross-section shape and the high-grade steel stress-strain 
relationship have a relatively minor effect on the inelastic 
behavior (Cohn, 1982, ref. 11). However, extensive 
experimental work is required in the future for a better 
evaluation of secondary moments in beams with a 
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fhe Design method 

A method for determining the ultimate load capacity of 
a continuous beam is proposed, which takes into account the 
secondary moment described by equation 5.1. Only a two-span 
symmetrical beam is considered, a similar approach being 
possibly applied to more complex structures. 

This method first evaluates the value of the net 
reinforcement index (w + wp -w') at both critical sections. 
The limiting value of 0.20 will be used as required by 
section 18.10.4 of the ACI 1977 Code to allow a variation of 
design moments from the elastic analysis. The limitation of 
0.30 is the dividing line between underreinforced and 
overreinforced members and this is confirmed by a parametric 
study (Cohn, 1982, ref. 11). Five different situations may 
arise as shown in Table 5.6. In cases 1, 2 and 3 the 
negative moments calculated by elastic analysis for any 
loading arrangement, may be increased or decreased by not 


more than 


W - Wp, - W 
20% SA tne DETCANE 
0.30 


In cases 4 and 5, no modification to the negative moment 
obtained from an elastic analysis is permitted. 

Cases 1 and 2: Since both critical sections are 
underreinforced, full redistribution is likely and hence no 
secondary moment remains at ultimate. The ultimate moments 


are those due to gravity loads only; 
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Moment used 
in 
design 


Modified 
moments 
according to 
Section 18.10.4 
of ACI 1977 


Elastic moments 
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Case 3: The section at maximum positive moment is 
overreinforced, which determines the ultimate load due to 
reduced ductility at that section. This ultimate load is 
calculated using the ratio of moment capacity provided, and 
then substituted into equation 5.1 to calculate the 
magnitude of secondary moment remaining at ultimate; 

Case 4: Although the beam is designed according to the 
elastic analysis, the low value of the reinforcement index 
at the support section allows the ultimate moments to be 
reached at all critical sections as in cases 1 and 2. This 
is achieved by yielding of the reinforcement and extensive 
cracking. The secondary moment is reduced to zero at 
ultimate; 

Case 5: The support section is overreinforced, and as a 
result the ultimate load corresponds to the load at which 
the support ultimate moment is attained through an elastic 
analysis. The procedure is similar to the one used in 
case 3. | 

The approach adopted in this design method satisfies at 
first the equilibrium condition by calculating the load from 
equations of statics, secondly the safety condition by 
ensuring that no premature failure occurs at any section of 
the beam, and thirdly the ductility condition to allow 
redistribution of forces from the elastic distribution to 
the assumed distribution. These three conditions of member 
proportioning correspond to the so-called lower bound 


solution, i.e. the structure is certainly capable to carry 
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the calculated load. 


As the member strength is approached the inelastic 
behavior at some sections results in a redistribution of 
moments. Recognition of this behavior can be advantageous, 
but a rigorous design method for moment redistribution is 
quite complex. The recommended design procedure is an 
attempt to give an account of the dependence of the 
secondary moment on the degree of cracking and moment 
redistribution. The amount of adjustment is kept within safe 


limits defined by current Code provisions. 
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6. SUMMARY, CONCLUSION AND RECOMMENDATIONS 


6.1 Summary 

The overall objective of this study was to investigate 
the flexural behavior characteristics of continuous 
prestressed concrete beams and determine the effects of 
inelastic behavior on the secondary moment. 

The analysis considered the equilibrium and strain 
compatibility of each beam section to evaluate the 
conditions at a particular cross-section of the beam. In 
addition, a theoretical stress-strain curve was used for 
both concrete and prestressing steel to establish the 
moment-average curvature relationships. Conventional 
principles of geometry were applied to determine the 
distribution of moment at any location in the beams. The 
complete moment-load curves to failure for several beams 
were obtained from the analysis and compared with available 


test results. 


6.2 Conclusion 

The reinforcement index value of 0.30 set by the ACI 
Building Code to define underreinforced sections can be used 
as a criterion to evaluate the ductility of continuous 
beams. For ductile beams, full redistribution of moment is 
possible due to yielding of the reinforcement and extensive 
cracking, with the understanding that minimum bonded 


reinforcement is provided to prevent the early failure of 
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the reinforcement. For non-ductile members, the ultimate 
limit state is defined by the load at which the moment 
capacity of a brittle section is attained. In such cases 
only partial moment redistribution is achieved. 

The secondary moment varies due to the readjustment of 
reactions caused by cracking, and this variation is 
dependent upon the redistribution of moment carried out at 
ultimate. A cubic variation is assumed in the present study, 
but a different type of variation may be set in the light of 


further investigation. 


6.3 Recommendation for future work 

More experimental programmes on continuous prestressed 
concrete beams are needed to study the phenomenon of 
secondary moment. A realistic evaluation of the secondary 
moment is necesSary because its effect on the load carrying 
Capacity of a given beam may be favorable or unfavorable, 
depending on the location of a particular cross-section. 

The analysis described in this investigation can be 


extended to more general statically indeterminate systems. 
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General Remarks 


This Appendix gives detailed explanation on the 


- computer program derived from the theory described in 


Chapter 4. The program consists of 8 subroutines and 


7 functions coordinated by the Main Routine. The sequence of 


execution is summarized in the diagram shown in the 


following page. The procedure of each subprogram is 


explained, and the program output for 6 beams is listed at 


the end of the Appendix. 


The limitations are: 
Two-Span symmetrical continuous prestressed beam with 
rotations allowed at all supports; 
Two symmetrical concentrated loads being applied at any 
position; 
Linear variation of tendon profile; 
Rectangular or doubly-symmetric I-shaped cross-section; 
No reinforcement other than prestressing bonded strands; 
The stress-strain relationships are to be modelized for 
both concrete and steel; only Grade 250 or 270 strand 
are considered; 

Secondary moments are neither evaluated nor included in 
the calculations; 


Imperial units. 
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MAIN ROUTINE 


The MAIN routine defines all loading stages starting 
from the load at which first cracking occurs, up to 
ultimate. The cross-section properties, the ultimate 
capacity and the cracking moment are evaluated at the 
critical sections by calling particular functions, after all 
the necessary informations are read. 

The elastic moment distribution along the beam can be 
calculated knowing the position of application of the 
concentrated load. If the moment capacity ratio provided at 
the critical sections is different from the elastic ratio, 
then redistribution of moment is likely. 

The MAIN routine assumes an initial set of moments and 
calls subroutine NEWTN so that this set of moments is varied 
until the compatibility equation is satisfied. The correct 
distribution of moments is then printed and another 
calculation is done for a higher load stage. 

At first cracking at a critical section, the beam is 
Still in the linear range and the compatibility 
equation 4.39 is used. All subsequent load stages are tested 


using compatibility equation 4.40. 
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NEWTN 


Newton-Raphson iterative procedure to evaluate a 
distribution of moments at critical sections that satisfies 
the compatibility conditions. The error margins adopted are: 

a. Elastic case, equation 4.39: + 0.0001 radians; 
b. Non-linear case, equation 4.40: + 0.05 radians. 

If x represents the distribution of moments and f(x) 
the deformations resulting from that distribution of 
moments, expressed in the compatibility equation, then the 
iterations are performed as shown in the following figure. 
The initial guess must be close enough to the solution in 


order to obtain convergence. 
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SEGM 


defines the portions of beam as shown in the following 
page. Beam segment lengths are chosen shorter in regions 
where the curvature reaches peak values. The distance to the 
end support, the eccentricity of steel and the moment at the 


center of each beam segment are evaluated. 


EACH 

evaluates the curvature and the contribution of each 
segment to the total deformations expressed in 
equations 4.39 and 4.40. If the moment at a segment is 
smaller than the cracking moment defined in equation 4.11, 
linear extrapolations are made to obtain the actual moment 
and curvature. Otherwise the general cracked section 
analysis of Section 4.3.4 is carried out to evaluate the 
steel stress at a crack. The average curvature in the 


segment in then calculated as indicated in Section 4.3.5. 
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CRITIM 


evaluates the cracking moment at a particular 


cross-section as described in Section 4.3.3. 


UMOM 

calculates the ultimate moment at a particular 
cross-section. If the net reinforcement index exceeds 0.30, 
a different equation is used as pointed out in 


Section 18.8.2 of the 1977 ACI Code. 


NLSYST 

Newton-Raphson iterative method with 2 parameters to 
solve 2 non-linear simultaneous equations. The derivatives 
of the equations considered have a contribution from each 
Parameter. The simultaneous equations are defined in 


Subroutine FCN. 


ELIM 
Gaussian elimination to calculate the corrections to be 


included in each iteration in Subroutine NLSYST. 
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FCN 


defines the equations of equilibrium and strain 


compatibility 4.5. 


AFCRAK 
General cracked section analysis defined in 
Section 4.3.4. The steel stress is calculated from 


equation 4.19. 


AXNT 
defines the equations of position of the neutral 


axace4..i2. and, 4.33. 
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PREL 
defines the stress-strain relationship for steel. The 
strain is expressed in terms of steel stress. 
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Stress-strain relationship for strands 


ROMB 

Romberg integration; performs the Simpson's rule for 
known values of a function. Intervals are then halved and 
results are extrapolated. The integration has been done on 
the curvatures obtained at 60 locations away from a crack to 


Calculate the average curvature defined in equation 4.35. 
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CMOM 
computes the contribution to the total moment, of the 
concrete area stressed in compression. The calculation takes 


into account the stress-strain relationship of concrete. 
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CMOM (cont'd) 


Numerical integration (Gaussian integration) 
The calculations in Functions CMOM and COMP are based 


on this method. The integration using 2 sample points 


located at t=-- 


exactly a cubic polynomial. 
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COMP 
computes the contribution to the total force, of the 


area stressed in compression. The procedure of integration 


is Similar to CMOM's. 
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COMP (cont'd) 
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PROGRAM TO ANALYZE THE MOMENT REDISTRIBUTION IN 
TWO-SPAN SYMMETRICALLY LOADED CONTINUOUS 
PARTIALLY PRESTRESSED CONCRETE BEAMS. 
CONCENTRATED SYMMETRICAL VERTICAL LOADS: 
STEEL TENDON PROFILE. 


LINEARLY VARYING 


MAIN PROGRAM 


DIMENSION X(2), F(2) 

COMMON X, F, B, H, EXL, EXS, PRESF, AREAP, CSTR, ELMOD, 

+FCU, E1, E2L, E2S, EPSO, SIGPU, AREAC, EPSU, COMPAT,DELTA(2), 
+ CTEN, ETEN, E2X, CRMOM, FS(60), STRAIN(6O), RESULT, 
+ FRACL, SPAN, ISECT, BW, HF, R2, FSB, CURVCT, ASR, ETOP, 
+ RDFACT, RATMOD, ALPHA, BETA, GAMMA, AXF, EXTMOM, TC1, 
+ ANGLEA, ANGLEB 


INPUT OF BEAM, CROSS-SECTION, AND MATERIAL DATA 
PARAMETERS ARE 


NUMBER REFERENCE NUMBER OF THE BEAM 

B OVERALL WIDTH OF CROSS-SECTION IN INCHES 

H OVERALL HEIGHT OF CROSS-SECTION IN INCHES 

EXL TENDON ECCENTRICITY AT LOAD POINT IN INCHES 

EXS TENDON ECCENTRICITY AT SUPPORT IN INCHES 

PRESF APPLIED PRESTRESSING FORCE IN LBS 

AREAP TOTAL AREA IN SQ. INCHES, OF PRESTRESSING TENDON 

SIGPU GRADE OF HIGH-STRENGTH PRESTRESSING STEEL IN KSI 

ELMOD MODULUS OF ELASTICITY IN PSI OF HIGH-STRENGTH STEEL 

CSTR CONCRETE STRENGTH IN PSI 

EPSO CONCRETE STRAIN AT MAXIMUM STRESS (MODEL) 

EPSU CONCRETE STRAIN AT CRUSHING (MODEL) 

RDFACT REDUCTION FACTOR OF STRENGTH OF CONCRETE (MODEL) 

SPAN SPAN LENGTH IN INCHES BETWEEN TWO SUPPORTS 

FRACL FRACTION OF SPAN FROM END SUPPORT WHERE ONE OF THE 
: CONCENTRATED LOADS IS APPLIED 

ISECT = O FOR A RECTANGULAR CROSS SECTION 

= 1 FOR AN I-SHAPED CROSS SECTION 
BW WEB THICKNESS IN INCHES (I-SECTION) 
HE FLANGE HEIGHT IN INCHES (I-SECTION) 


qgagqagqaqagqNaganaanaNnananaNnaanaNnanaNanannanananna 


READ(5, 1000) 


1000 FORMAT (110) 
READ(5, 2000) 

2000 FORMAT (4F10. 
READ(5, 2001) 

2001 FORMAT (3F10. 
-  READ(5, 2002) 
2002 FORMAT (4F10. 
READ(5, 2003) 

2003 FORMAT(2F10.5, 


IF(ISECT .EQ. 
READ(5, 2004) 


NUMBER 
B, H, EXS 
5) 

PRESF, AREAP, 
5, E10.5) 
EPSO, 


EXL, 


SIGPU, ELMOD 


CSTR, EPSU, RDFACT 
5) 
SPAN, 

110) 
0) GO TO 10 
BW, HF 


FRACL, ISECT 


FORMAT(2F10.5) 
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Cc 
C ECHO CHECK OF INPUT DATA 
C 
POS = SPAN*FRACL 
WRITE(6,500) NUMBER,POS 
IF (ISECT .EQ. O) WRITE(6,501) B, H 
WRITE(6,502) B, H, BW, HF 
WRITE(6,600) EXL, EXS, SPAN 
WRITE(6,700) CSTR, EPSO, EPSU, RDFACT 
WRITE(6,800) AREAP, PRESF, SIGPU, ELMOD 
500 FORMAT(50(’*’),/’MOMENT-LOAD CURVE FOR THE TWO-SPAN’, 
+ ‘ CONTINUOUS PRESTRESSED BEAM NUMBER ’,I5,’ LOADED’, 
+ /’SYMMETRICALLY WITH TWO CONCENTRATED VERTICAL LOADS’, 
+ F7.2,’ INCHES FROM THE END SUPPORTS. ’) 


501 FORMAT(/’RECTANGULAR CROSS-SECTION, WIDTH = ‘’,F7.2, 
+ 4+IN«;, HEIGTH = ’,F7.2,’ IN.’) 

502 FORMAT(/’I-SHAPED SECTION, WIDTH = ’,F7.2, 
+ ’ IN@epgyHELGHT =e’ sk7e2,4eENa¢, 
+ P2Q2X ee WEG rd ce Ns 5 FLANGE 2) F7.2,7 JIN. 2) 


600 FORMAT(//’‘TENDON ECCENTRICITY = ’,F10.2,’ IN. AT LOAD POINT’ 


+ wif22%, £1062, ~ BN. AT SUPPORT. ’, 
+ /’SPAN BETWEEN SUPPORTS ’,F10.2,’ IN.’) 
700 FORMAT(//’CONCRETE STRENGTH ’,F15.2,’ PSI ,’, 
¥ fe STRAIN EPSO ’,F15.5, 
+ /' EPSU ’,F15.5, 
+ , REDUCTION FACTOR ’,F10.2) 
800 FORMAT(//’TOTAL TENDON AREA ’,F15.4,’ SQ. IN.’, 


+ /’APPLIED PRESTRESSING’ ,F13.2,’ LBS’, 
+ /’GRADE False en AK SS, 
ot /’MOD. OF ELASTICITY ’,€1939.3,’ PSI’, 
+ /50(’*’),//’THE CRACKING MOMENTS AT SUPPORT AND AT THE’, 
+ ’ LOAD POINT ARE, RESPECTIVELY: ’) 
Cc 
C CALCULATIONS OF SECTION PROPERTIES 
Cc 
10 FCU = CSTR/(.8 + .0001*CSTR) 
RATMOD = ELMOD/(57000. *SOQRT(CSTR) ) 
ALPHA = (B-BW)*HF + RATMOD*AREAP 
E1 = PRESF/(AREAP*ELMOD ) 
C 
Cc FCU AVERAGE CONCRETE STRESS (FROM WARWARUK, REF. 28) 
C RATMOD RATIO OF STEEL MODULUS / CONCRETE MODULUS 
C ALPHA COEFFICIENT TO BE USED IN FUNCTION "AXNT" 
Cy E INITIAL STEEL STRAIN (STAGE 1) 
C 
IF (ISECT .EQ. O) GO TO 12 
I MOMENT OF INERTIA OF CROSS-SECTION 
AREAC AREA OF CONCRETE 
R2 SQUARED RADIUS OF GYRATION 


GYRL,GYRS TERMS FROM EQ. 3.3 
E2L,E2S STEEL STRAIN AT STAGE 2 
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I = 2.*(B*HF**3/12.) + .5*B*HF*(H-HF) **2 


+ + BW*(H-2. *HF)**3/12. 
AREAC = 2.*B*HF + BW*(H-2.*HF) 
R2 = I/AREAC 
GYRL = EXL*EXL/R2 
GYRS = EXS*EXS/R2 
GO TO 13 
12 AREAC = B*H 
GYRL = 12.*(EXL/H)*(EXL/H) 
GYRS = 12.*(EXS/H)*(EXS/H) 
13 E2L = PRESF*(1+GYRL)/(AREAC*57000. *SQRT(CSTR) ) 
E2S = PRESF*(1+GYRS)/(AREAC*57000. *SQRT(CSTR) ) 


CALCULATION OF ULTIMATE MOMENTS (ITERATIVE PROCEDURE ) 


X(1) 
X(2) 


X(1) 
xX(2) 
UMOML 
UCURVL 


X(1) = 
X62), o* 
UMOMS 
UCURVS 


CALCULATI 
CTEN 
ETEN 


ETEN = 


EVALUATIO 
ELAS 
PROVID 


EL 
PR 
IF 
IF 


EVALUATIO 


50 WR 
CA 
GO 
60 WR 
CA 


FIRST GUESS OF LOCATION OF NEUTRAL AXIS 
FIRST GUESS OF STEEL STRESS 


(EXL+H/2.)/10. 

SIGPU- 30. 
= UMOM(1, EXL) 

= (PREL(X(2)) + EPSU)/(EXL+H/2. ) 


(EXS+H/2.)/10. 

SIGPU - 30. 
= UMOM(2,EXS) 

= (PREL(X(2)) + EPSU)/(EXS+H/2. ) 


ON OF CRACKING MOMENTS 
MODULUS OF RUPTURE OF CONCRETE 
CONCRETE STRAIN AT TENSILE STRENGTH 


7.5*SQRT(CSTR) 
7.5/57000. 

CRITIM(EXS) 
CRITIM(EXL) 


N OF MOMENT RATIOS 
ELASTIC MOMENT RATIO 
MOMENT CAPACITY RATIO PROVIDED 


AS = (1.-FRACL)*(2.+FRACL)/(1.+FRACL) 
OVID = UMOML/UMOMS 
(ABS(ELAS-PROVID) .LE. 0.1) GO TO 80 
(ELAS*CSMOM - CLMOM) 50, 50, 60 


N OF THE FIRST CRACKING CONDITIONS 


ITE (6,6000) 

LL NEWTN (0, ELAS*CSMOM, CSMOM, 1) 
TO 70 

ITE (6,6001) 

LL NEWTN (0, CLMOM, CLMOM/ELAS, 1) 
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LOADING STAGES BETWEEN FIRST CRACKING AND ULTIMATE 


70 DO 51 KS = 1,5 
ADD = 0.0 
DO 51 KSS = 1,5 
SM = KS*(1.+ADD)*CSMOM 
PM = PROVID*SM 
IF ((SM.GT.UMOMS) .OR. (PM.GT.UMOML)) GO TO 30 
WRITE(6,4000) SM, PM 
ADD = ADD + 0.2 
51 CALL NEWTN (0, PM, SM, 2) 


ULTIMATE STAGE 


30 WRITE(6,3000) UMOMS, UCURVS, UMOML, UCURVL 
CALL NEWTN (0, UMOML, UMOMS, 2) 
STOP 
4000 FORMAT(/50(’-’),/’ THE NEXT ASSUMED SET OF MOMENTS IS:’, 
+ /15X,’‘AT SUPPORT : ’,£13.6,’ ,AND AT LD POINT : “, 
+UETS.067°" CE- IN’) 
3000 FORMAT (/,50(’*’),/’ULTIMATE STAGE: ASSUMED’, 


+ /’ SUPPORT MOMENT = ’.£13.6,’ LB-IN.’, 
+ i: CURVATURE = ',£13.6, 

6 /’ LOAD POINT MOMENT =’,£13.6,’ LB-IN.’, 
+ , CURVATURE = ’,£13.6) 


6000 FORMAT(/,50(’-’),/’CRACKING OCCURS FIRST AT SUPPORT’ ) 
6001 FORMAT(/,50(’-’),/’CRACKING OCCURS FIRST AT LOAD POINT’) 


STOP 
IF THE MOMENT CAPACITY RATIO PROVIDED IS EQUAL TO THE ELASTIC 
ONE, THEN NO REDISTRIBUTION OF MOMENT OCCURS. PRINT MESSAGE. 
80 WC = 2.*(CLMOM+CSMOM*FRACL )/(FRACL*(1.-FRACL)*SPAN) 
WU = 2.*(UMOML+UMOMS*FRACL )/(FRACL*(1.-FRACL) *SPAN) 
WRITE (6,8000) PROVID,ELAS, WC,CSMOM,CLMOM, WU,UMOMS , UMOML 
8000 FORMAT (/’CAPACITY PROVIDED, ’,F7.2,’ ACCORDING TO’, 
=f ‘ ELASTIC RATIO ’,F7.2,’ :* NO REDISTRIBUTION. *’, 
+ //'CRACKING : OCCURS AT LOAD = ’,£13.6,’ LBS’, 
+” Lé SUPPORT MOMENT = ’,E13.6,’ LB-IN.’, 
+ ff LD PT MOMENT eee EG. Orem bal N ee, 
+ //‘ULTIMATE : OCCURS AT LOAD = ’,E13.6,’ LBS’, 
+ he SUPPORT MOMENT = ’,E13.6,’ LB-IN.’, 
+ /i LD PT MOMENT ee fio. 6 el RINE) 
STOP 
END 
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SUBROUTINE NEWTN (I, PLMOM, SPMOM, KANG) 
SUBROUTINE THAT PERFORMS THE NEWTON-RAPHSON ITERATIVE 
METHOD TO VARY EITHER THE LOAD POINT MOMENT OR THE 
SUPPORT MOMENT UNTIL SATISFACTION OF THE COMPATIBILITY 
EQUATION. 
PARAMETERS ARE - 

I O IF EACH INTERMEDIATE ITERATION IS TO BE PRINTED 
1 IF NOT 


PLMOM NITIAL ASSUMED MOMENT AT LOAD POINT 
SPMOM NITIAL ASSUMED MOMENT AT SUPPORT 
KANG 1 INDICATES ELASTIC COMPATIBILITY EQUATION 


“nea HW ot 


2 INDICATES NON-LINEAR COMPATIBILITY EQUATION 


DIMENSION X(2), F(2) 

COMMON X, F, B, H, EXL, EXS, PRESF, AREAP, CSTR, ELMOD, 

+FCU, E1, E2L, E2S, EPSO, SIGPU, AREAC, EPSU, COMPAT,DELTA(2), 

+ CTEN, ETEN, E2X, CRMOM, FS(60), STRAIN(60), RESULT, 

FRACL, SPAN, ISECT, BW, HF, R2, FSB, CURVCT, ASR, ETOP, 
RDFACT, RATMOD, ALPHA, BETA, GAMMA, AXF, EXTMOM, TCi, 
ANGLEA, ANGLEB 


t+ + 


SET TOLERANCE VALUES FOR X, F(X) TO TERMINATE ITERATIONS. 


NLIM : LIMIT TO NUMBER OF ITERATIONS 
XTOL = 9000. 

FTOL = .0001 

IF (KANG,.EQ. 2) FTOL. = 3. 

NLIM = 10 


DELTM = 10000. 

LOGICAL PRINT 

PRINT = .TRUE. 

IF (I .NE. O) PRINT = .FALSE. 
CALL SEGM (PLMOM, SPMOM, KANG) 
FX = COMPAT 


ITERATIONS ARE PERFORMED ON THE SUPPORT MOMENT FIRST 


CALL SEGM (PLMOM, SPMOM+DELTM, KANG) 
FXD = COMPAT 


FDER = (FXD - FX)/DELTM 
DELX = FX/FDER 
IF (ABS(DELX) .GE. 10000.) DELX = -10000. 


SPMOM2 = SPMOM - DELX 
SPMOM1 = SPMOM 
CALL SEGM (PLMOM, SPMOM2, KANG) 
FX = COMPAT 
KU = O 
IF (.NOT. PRINT) GO TO 9 
WRITE (6,199) JU, KU, PLMOM, SPMOM2, FX 
9 IF (KJ .EQ. 0) SPMOM = SPMOM2 
IF (ABS(DELX) .LE. XTOL) GO TO 60 
IF (ABS(FX) .LE. FTOL) GO TO 70 
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IF TOLERANCE IS NOT MEET, ITERATIONS ARE DONE ON THE 
LOAD POINT MOMENT. 


SPMOM = SPMOM1+ 1000. 
DO 20 KU = 1, 8 
PLMOMD = PLMOM + DELTM 
CALL SEGM (PLMOMD, SPMOM, KANG) 
FXD = COMPAT 
FDER = (FXD - FX)/DELTM 
DELX = FX/FDER 
IF (ABS(DELX) .GE. 10000.) DELX = -10000. 
PLMOM = PLMOM - DELX 
CALL SEGM (PLMOM, SPMOM, KANG) 
FX = COMPAT 
80 IF (.NOT. PRINT) GO TO 5 
WRITE (6,199) JU, KJ, PLMOM, SPMOM, FX 
5 IF (ABS(DELX) .LE. XTOL) GO TO 60 
IF (ABS(FX) .LE. FTOL) GO TO 70 
20 CONTINUE 


WHEN LOOP IS NORMALLY COMPLETED, NLIM IS EXCEEDED. 


WRITE (6,200) NLIM, PLMOM, SPMOM, FX 
RETURN 


THIS SECTION RETURNS AFTER MEETING TOLERANCE ON XTOL 
(MOMENT DISTRIBUTION IS SATISFACTORY) 


60 WRITE (6,202) JU, KJ, PLMOM, SPMOM, FX 
GO TO 300 


THIS SECTION RETURNS AFTER MEETING F(X) TOLERANCE 
(COMPATIBILITY EQUATION IS SATISFIED) 


70 WRITE (6,203) JU, KJ, PLMOM, SPMOM, FX 


CALCULATION OF LOAD CORRESPONDING TO FINAL DISTRIBUTION 
OF MOMENTS 


300 W = 2.*(PLMOM+SPMOM*FRACL )/(FRACL*BB*SPAN) 


WRITE(6,205) W 


205 FORMAT(//5X,’THE TOTAL APPLIED LOAD IS= ’,E13.6,’ LBS’) 


199 FORMAT (/5X,’AT ITER. ’,14,’ ’,14,’ THE MOMENTS AT ’ 
+ ‘LD POINT = ’,£13.6,’ AND AT SUPPORT = ’,£13.6, 
+ ’ LB-IN.’,/30X,’*** COMPATIBILITY EQUATION = ’, E13.6) 
200 FORMAT (//5X,’TOLERANCE NOT MET. AFTER ’, 14,’ ITER.’, 
+ ‘ LD PT. MOM. = ’, E13.6,’ SUP = ’, E13.6, 
+ ’ LB-IN.’,/30X,’*** COMPATIBILITY EQ. = ’,£13.6) 
202 FORMAT (//5X,’MOM. TOLERANCE MET IN ’, 14,’ ’,14,’ ITER.’, 
eet LOPES MOMS =) EIS; Goo VSUPL ear Eda. a 
+ ’ LB-IN.’,/30X’*** COMPATIBILITY EQ. = ‘,£13.6) 
203 FORMAT (//5X,’COMPAT. TOLERANCE MET IN ’, 14,’ ’,14, 
i @ETTERS. “LD PE.) MOM: =o *306819.6)77. .,.5UP. = 5 elo. 
+ *’ LB-IN.’,/30X,’*** COMPATIBILITY EQ. = ’,&13.6) 
RETURN 
END 


101 


rot 


Gate iu ba 
ee 


on 
eres a @ 


i 
a 


SH) WO aM00 3RA 1 sete 


O01 2 Kae, ‘A. af 3 
me Fuso - meses 
Lue ' Tr MONI9) MORE 
a) taemer « 
woke 2 ‘OD 1WALAT -YOu.) FT 
mae boty um & fh .9) 2tis 


a ly 
. mh 
Sabb OMe th 


3359¥9 @i mat cae | J 4 SOO 21 a00 ¥ 
A mom, tie a. czech" vm ce 
in sigciaae 
1oTx Ho Toe IIGT pvr ' 93M garne a ates vorrade rate 
(VEST IAAL GAR BY" Novem Ara To Taso) 


£1 eer one ad * a 3 ao a 


Sova a uOm, Ox 4 bedi is a Bore a 3 
vi! na seas 3 


i 


j 


GST: 1SeT Ae oy 


: if K2 MOMEE MOMS Got LU nmr | 
VO! PUTRT SIO! WARS of aapretaeee nie uA pian > 
: r a : Pp al 7? ™ = 
a a ae : | » (MASBe age Pep aaiee: or Air a ae 


(‘25% ‘28 aoa ee? Gaoo 03%, sna + 
“, 7 h- 2) ers at 4 ate 
: ; a.eb3ip a? AAG 

(aera HOST aus yt1408 


“ far “igs 


oy 


2 MS tee 4 ie th sf 
“or 17 afr 


aqaan0n 


qgaagagagaaa 


gagagananAa 


aaqgqagaa 


SUBROUTINE SEGM (PLMOM, 


SPMOM, KANG) 
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SUBROUTINE DISCRETIZING THE BEAM INTO SEGMENTS AS SHOWN IN 


PAGE 82. 
TENDON ECCENTRICITY, 


EVALUATED. 
THE SUBROUTINE EACH IS THEN CALLED. 
PARAMETERS ARE IDENTICAL TO THE SUBROUTINE NEWTN’S. 


LOCATION WITH RESPECT TO THE END SUPPORT, 
VALUE OF MOMENT OF EACH SEGMENT IS 


DIMENSION X(2), F(2), 2ZL(6), ZR(6), ZRS(6) 

COMMON X, F, B, H, EXL, EXS, PRESF, AREAP, CSTR, ELMOD, 

+FCU, E1, E2L, E2S, EPSO, SIGPU, AREAC, EPSU, COMPAT,DELTA(2), 
+ CTEN, ETEN, E2X, CRMOM, FS(60), STRAIN(60), RESULT, 

+ FRACL, SPAN, ISECT, BW, HF, R2, FSB, CURVCT, ASR, ETOP, 

ct: RDFACT, RATMOD, ALPHA, BETA, GAMMA, AXF, EXTMOM, TCi, 

+ ANGLEA, ANGLEB 


INITIALIZATION OF CUMULATIVE VARIABLES. 


COMPAT VALUE OF COMPATIBILITY EQUATION 
ANGLEA ANGLE AT END SUPPORT 
ANGLEB ANGLE AT ONE SIDE OF CENTRE SUPPORT 
COMPAT = O. 
ANGLEA = O. 
ANGLEB = O. 


PORTION OF BEAM BETWEEN END SUPPORT AND LOAD APPLICATION POINT. 
DISTANCE OF THE SEGMENT CENTER TO THE END SUPPORT 


200 


DISTM 
SLENGF 


ZL(1) 
ZL(2) 
ZL(3) 
ZL (4) 
ZL(5) 
ZL(6) 


ECCENTRICITY OF STEEL TENDON 
VALUE OF MOMENT 
SEGMENT LENGTH 


FRACL*SPAN/4. 
2.2*ZL(1) 
2.6*ZL(1) 
3.0*ZL(1) 
3.4*ZL(1) 
3.8*ZL(1) 

DO 200 IDIST=1,6 
EZ = EXL*ZL(IDIST)/(4.*ZL(1)) 
DISTM = PLMOM*ZL(IDIST)/(4.*ZL(1)) 
SLENGF = 2.*ZL(1) 

IF (IDIST .EQ. 1) GO TO 200 
SLENGF = 0.4*ZL(1) 
CALL EACH(ZL(IDIST), 


EZ, DISTM, SLENGF, 


KANG) 
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~ PORTION OF BEAM BETWEEN LOAD APPLICATION POINT AND THE POINT 
_WHERE THE TENDON ECCENTRICITY MEETS THE CROSS SECTION 
CENTER OF GRAVITY. 


300 


VAL 
ZR(1) 
ZR(2) 
ZR(3) 
ZR(4) 
ZR(5) 
ZR(6) 


DISTANCE OF SEGMENT TO LOAD APPLICATION POINT 
TENDON ECCENTRICITY 


DISTM | VALUE OF MOMENT IN SEGMENT 


(1-FRACL ) *SPAN 

VAL*EXL/(20.*(EXL+EXS) ) 

3.*ZR(1) 

5.*ZR(1) 

7.*ZR(1) 

9.*ZR(1) 

15.¥*ZR(1) 

DO 300 IDIST=1,6 

ER = EXL - EXL*ZR(IDIST)/(20. *ZR(1)) 

DISTM = PLMOM - ZR(IDIST)*(PLMOM+SPMOM) /VAL 
SLENGF = (2./3.)*ZR(6) 

IF (IDIST .EQ. 6) GO TO 300 

SLENGF = ZR(6)/7.5 

CALL EACH(4.*ZL(1)+ZR(IDIST), ER, DISTM, SLENGF, KANG) 


REMAINING PORTION OF BEAM (UP TO CENTRE SUPPORT) 


400 


DISTANCE OF SEGMENT TO THE LOAD APPLICATION POINT 
TENDON ECCENTRICITY 


ZRS(1) = VAL*EXS/(4.*(EXL+EXS)) 
ZRS(2) = 2.2*ZRS(1) 
ZRS(3) = 2.6*ZRS(1). 
ZRS(4) = 3.0*ZRS(1) 
ZRS(5) = 3.4*ZRS(1) 
ZRS(6) = 3.8*ZRS(1) 
PINF = EXL*VAL/(EXL+EXS) 
DO 400 IDIST=1,6 
ERS = -EXS*ZRS(IDIST)/(4.*ZRS(1)) 


DISTM = PLMOM - (PINF + ZRS(IDIST))*(PLMOM + SPMOM)/VAL 
SLENGF = 2.*ZRS(1) 

IF (IDIST .EQ. 1) GO TO 400 

SLENGF = 0.4*ZRS(1) 

CALL EACH(4.*ZL(1)+PINF+ZRS(IDIST), ERS,DISTM, SLENGF ,KANG) 


THIS SECTION EVALUATES THE COMPATIBILITY EQUATION AFTER 
THE CONTRIBUTION OF ALL THE SEGMENTS TO THE TOTAL ANGLES 
HAS BEEN ACCUMULATED. 


IF (KANG.EQ.2) COMPAT = (180./ 3.141593) * (FRACL. * 


(ABS(ANGLEA) + ABS(ANGLEB)) - ABS(ANGLEB) ) 


RETURN 
END 
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SUBROUTINE EACH(Z, EZ, DISTM, SLENGF, KANG) 


SUBROUTINE FOR CALCULATION OF CONTRIBUTION OF EACH SEGMENT 
CURVATURE TO THE TOTAL ANGLES IN THE COMPATIBILITY EQUATIONS. 
PARAMETERS ARE - 

Z DISTANCE OF CENTRE OF SEGMENT TO THE END SUPPORT 

EZ ECCENTRICITY OF SEGMENT 

DISTM VALUE OF MOMENT 

SLENGF LENGTH OF SEGMENT 

KANG =1 FOR ELASTIC COMPATIBILITY EQUATION 

=2 FOR NON-LINEAR COMPATIBILITY EQUATION 


DIMENSION X(2), F(2) 

COMMON X, F, B, H, EXL, EXS, PRESF, AREAP, CSTR, ELMOD, 

+FCU, E1, E2L, E2S, EPSO, SIGPU, AREAC, EPSU, COMPAT,DELTA(2), 
+ CTEN, ETEN, E2X, CRMOM, FS(60), STRAIN(60), RESULT, 

+ FRACL, SPAN, ISECT, BW, HF, R2, FSB, CURVCT, ASR, ETOP, 

+ RDFACT, RATMOD, ALPHA, BETA, GAMMA, AXF, EXTMOM, TCi, 

+ ANGLEA, ANGLEB 


IF THE MOMENT IN SEGMENT IS SMALLER THAN THE CRACKING 
MOMENT CORRESPONDING TO THE ECCENTRICITY, THEN ITS 
ACTUAL MOMENT IS EQUAL TO A PORTION OF THE CRACKING MOMENT. 


CRMOM = (DISTM/ABS(DISTM) )*CRITIM(EZ) 
IF (ABS(DISTM) .LT. ABS(CRMOM)) GO TO 5000 


GENERAL CRACKED SECTION ANALYSIS AS DESCRIBED IN 

CHAPTER 3. EVALUATION OF THE COEFFICIENTS TO BE USED 

IN THE EQUATION OF NEUTRAL AXIS. 

IF THE SEGMENT SECTION IS CRACKED, THE PRESTRESSING 

STEEL STRESS IS EVALUATED USING THE LINEAR BOND VARIATION 
CONCEPT. THE AVERAGE CURVATURE IS CALCULATED BY INTEGRATING 
THE CURVATURES THROUGH SUBROUTINE ROMB. 


BETA = (B-BW)*HF*HF/2. + RATMOD*AREAP*(ABS(EZ)+H/2. ) 
GAMMA = (B-BW)*HF**3/3. + RATMOD*AREAP*(ABS(EZ)+H/2.)**2 
AXF = AREAP*ELMOD*(E1+E2X) 


CONDITION AT THE CRACK, CRACKING MOMENT HAS BEEN JUST REACHED 


EXTMOM = ABS(CRMOM) - AXF*(ABS(EZ)+H/2. ) 
Ye=,H/2% 
CALL AFCRAK(Y, EZ, 1, LENGTH) 

FSA = X(2)/1000. 


CONDITION AT THE CRACK, CRACKING MOMENT EXCEEDED 


EXTMOM = ABS(DISTM) - AXF*(ABS(EZ)+H/2.) 
Ye2 H/2: 

CALL AFCRAK(Y, EZ, 1, LENGTH) 

FSCR = X(2)/1000. 
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ETOP IS THE DISTANCE FROM THE TOP OF THE SECTION TO THE 
"FICTITIOUS" FORCE 


ETOP = (EXTMOM - AXF*ABS(EZ))/AXF - H/2. 


THE PRESTRESSING STEEL STRESS VARIATION IS EVALUATED 
AT 60 POINTS EVENLY SPACED WITHIN 2 HYPOTHETICAL 
CRACKS IN THE CONCRETE. 


DLB = O. 
DO 80 LENGTH=1,60 
DLB = DLB + 0.0125 
FS(LENGTH) = FSCR -2.*(FSA-FSB)*(DLB-.5*DLB*DLB) 
Yeertifa:. 
CALL AFCRAK(Y, EZ, 2, LENGTH) 
STRAIN(LENGTH) = X(1) + PREL(FS(LENGTH) ) 
80 CONTINUE 
CALL ROMB(60, EZ, DISTM) 
GO TO 100 


5000 RESULT = (DISTM/ABS(CRMOM) ) *CURVCT 


THE AVERAGE CURVATURE FROM THE CRACKED SEGMENTS, 
OR THE CURVATURE FROM THE UNCRACKED SEGMENTS 
IS ASSIGNED TO THE VARIABLE "RESULT" 


100 IF (KANG .NE.2) COMPAT = COMPAT + RESULT*Z/SPAN 
ANGLEA = ANGLEA + RESULT*(SPAN-Z)/SPAN 
ANGLEB = ANGLEB + RESULT*Z/SPAN 
RETURN 
END 
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FUNCTION CRITIM (EZ) 
FUNCTION FOR THE EVALUATION OF THE CRACKING MOMENT OF 
EACH SEGMENT. THE SINGLE PARAMETER IS THE TENDON 
ECCENTRICITY IN THE SECTION. 
DIMENSION X(2), F(2) 
COMMON X, F, B, H, EXL, EXS, PRESF, AREAP, CSTR, ELMOD, 
+FCU, E1, E2L, E2S, EPSO, SIGPU, AREAC, EPSU, COMPAT,DELTA(2), 
+ CTEN, ETEN, E2X, CRMOM, FS(60), STRAIN(60), RESULT, 
+ FRACL, SPAN, ISECT, BW, HF, R2, FSB, CURVCT, ASR, ETOP, 
+ RDFACT, RATMOD, ALPHA, BETA, GAMMA, AXF, EXTMOM, TCi, 
+ ANGLEA, ANGLEB 
SET THE X INITIAL GUESSES 
X(1) = 0.0005 
X(2) = 100. 
SECTION PROPERTIES 
GYRZ : TERM OF EQUATION 3.3 
E2X : STEEL STRAIN AT STAGE 2 
IF (ISECT .EQ. 0) GO TO 21 
GYRZ = ABS(EZ)**2/R2 
GO TO 23 
21 GYRZ = 12.*(ABS(EZ)/H)*(ABS(EZ)/H) 
23 E2X = EPSO*(1-SQRT(1-PRESF/(AREAC*CSTR) )) *(1+GYRZ) 
SIMULTANEOUS EQUATIONS OF EQUILIBRIUM AND STRAIN 
COMPATIBILITY ARE SOLVED THROUGH SUBROUTINE "NLSYST". 
CASE OF THE CRACKING CONDITIONS. 
CALL NLSYST (2, 1, 3, EZ, LENGTH) 
RESULTING STEEL STRESS AND NEUTRAL AXIS 
FSB = X(2) 
C = H*X(1)/(ETEN+X{(1)) 
IF(ISECT .EQ. 0) GO TO 24 
IF(H-C .LE. HF) GO TO 24 
CRITIM = ABS(CMOM(C,X(1))) + BW*CTEN*(H-C)*(H-C)/3. 
a + (B-BW) *CTEN*(H-C-HF/2. ) *HF 
+ + AREAP*X(2)*1000. *(ABS(EZ)+H/2.-C) 
GO TO 25 
24 CRITIM = ABS(CMOM(C,X(1))) + B*CTEN*(H-C)*(H-C)/3. + 
+ AREAP*X(2)*1000. *(ABS(EZ)+H/2.-C) 
25 CURVCT = (X(1)+ETEN)/H 
IF ((ABS(EZ) .EQ@. EXL) .OR. (ABS(EZ) .EQ. EXS)) GO TO 500 
RETURN 
500 WRITE(6,3000) CRITIM, CURVCT 
3000 FORMAT (/’ MOMENT =n? €49.6>* LB-IN 
+ /’ CURVATURE = ',E13.6) 
RETURN 
END 
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FUNCTION UMOM(MCASE, ECC) 


FUNCTION FOR CALCULATING THE ULTIMATE MOMENT. 
PARAMETERS ARE - 
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MCASE = 1 FOR THE EVALUATION OF LOAD POINT MOMENT 


= 2 FOR THE EVALUATION OF SUPPORT MOMENT 


ECC ECCENTRICITY OF PRESTRESSING STEEL TENDON 


DIMENSION X(2), F(2) 
COMMON xX, F, B, H, EXL, EXS, PRESF, AREAP, CSTR, 


ELMOD, 


+FCU, E1, E2L, E2S, EPSO, SIGPU, AREAC, EPSU, COMPAT,DELTA(2), 


+ CTEN, ETEN, E2X, CRMOM, FS(60), STRAIN(60), 


ttt 


ANGLEA, ANGLEB 


VALUES SET FOR USE IN SUBROUTINE "NLSYST" 


N = 2 

DELTA(1) = 0.001 
DELTA(2) = 2. 

EZa air 

LENGTH = 1 


CALL NLSYST (N, 1, MCASE, EZ, LENGTH) 
CHECK THE VALUE OF NET REINFORCEMENT INDEX 


OMEGA = AREAP*X(2)*1000./(B*(ECC+H/2.)*CSTR) 


TF) Ot FS. (OS) GOFTORTS 
IF (ISECT .EQ. 0) GO TO 40 
IF (X(2)*1000. *AREAP/(FCU*B) .LE. HF) GO TO 40 


FRACL, SPAN, ISECT, BW, HF, R2, FSB, CURVCT, ASR, 
RDOFACT, RATMOD, ALPHA, BETA, GAMMA, AXF, EXTMOM, TC1, 


RESULT, 
ETOP, 


DIFFERENT FORMULAS ARE TO BE USED WHEN THE REINFORCEMENT 


INDEX EXCEEDS 0.3 


IF (OMEGA - 0.30) 301,301,300 
300 WRITE (6,4000) OMEGA 
UMOM = (CSTR*BW*(ECC+H/2)**2)/4. 
+ + 0. 85*CSTR*(B-BW) *HF *(ECC+H/2.-HF/2.) 
RETURN 
301 ASR = AREAP - FCU*(B-BW) *HF/(X(2)*1000. ) 
UMOM = (AREAP-ASR) *X(2) * 1000. *(ECC+H/2.- .42*X(1)) 
+ + ASR*X(2)* 1000. *(ECC+H/2.-HF/2.) 
RETURN 
40 IF (OMEGA - 0.30) 201,201,200 
200 WRITE (6,4000) OMEGA 
UMOM =(CSTR*B*(ECC+H/2.)**2)/4. 
RETURN 
201 UMOM = AREAP*X(2)* 1000. *(ECC+H/2.-.42*X(1)) 
RETURN 
70 WRITE (6,7000) X(1) 


4000 FORMAT (/’OVERREINFORCED SECTION ***’, 


. /'****OMEGA = ’,£13.5,’ GT. oO. 207) 


STOP 
END 


7000 FORMAT (/’ NEGATIVE N. A. = £,£12:5,° CULTIM: STAGE)’) 
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SUBROUTINE NLSYST (N, I, MCASE, EZ, LENGTH) 


THIS SUBROUTINE SOLVES A SYSTEM OF N NON-LINEAR EQUATIONS BY 
NEWTON’S METHOD. THE PARTIAL DERIVATIVES OF THE FUNCTIONS ARE 
ESTIMATED BY DIFFERENCE QUOTIENTS WHEN A VARIABLE IS PERTUBED 
BY AN AMOUNT EQUAL TO DELTA (DELTA IS ADDED). 

THIS IS DONE FOR EACH VARIABLE IN EACH FUNCTION. INCREMENTS 
TO IMPROVE THE ESTIMATES FOR THE X-VALUES ARE COMPUTED FROM 

A SYSTEM OF EQUATIONS USING SUBROUTINE ELIM. 

PARAMETERS ARE - 


FCN SUBROUTINE THAT COMPUTES VALUES OF THE FUNCTIONS. 

N THE NUMBER OF EQUATIONS 

MAXIT LIMIT TO THE NUMBER OF ITERATIONS THAT WILL BE USED 

X ARRAY TO HOLD THE X VALUES. INITIALLY THIS ARRAY 
HOLDS THE INITIAL GUESSES. IT RETURNS THE FINAL 
VALUES. 

F AN ARRAY THAT HOLDS VALUES OF THE FUNCTIONS 

DELTA A SMALL VALUE USED TO PERTURB THE X VALUES SO 
PARTIAL DERIVATIVES CAN BE COMPUTED BY DIFFERENCE 
QUOTIENT 

XTOL TOLERANCE VALUE FOR CHANGE IN X VALUES TO STOP 


ITERATIONS. WHEN THE LATEST CHANGE IN ANY X 
MEETS XTOL, THE SUBROUTINE TERMINATES. 


FTOL TOLERANCE VALUE ON F TO TERMINATE. WHEN THE 
LATEST F VALUE IS LESS THAN FTOL, SUBROUTINE 
TERMINATES. 


DIMENSION X(2), F(2), A(2,3), XSAVE(2),FSAVE(2) 
COMMON X, F, B, H, EXL, EXS, PRESF, AREAP, CSTR, ELMOD, 
+FCU, E1, E2L, E2S, EPSO, SIGPU, AREAC, EPSU, COMPAT,DELTA(2), 
+ CTEN, ETEN, E2X, CRMOM, FS(60), STRAIN(60), RESULT, 
+ FRACL, SPAN, ISECT, BW, HF, R2, FSB, CURVCT, ASR, ETOP, 
“t- RDFACT, RATMOD, ALPHA, BETA, GAMMA, AXF, EXTMOM, TC1, 
+ ANGLEA, ANGLEB 


FTOL = 0.0005 
LOGICAL PRINT 
CHECK VALIDITY OF VALUE OF N 


IF (N.LT.2 .OR. N.GT.3) GO TO 999 
PRINT = .TRUE. 


I INDICATES IF INTERMEDIATE RESULTS ARE TO BE PRINTED 


IF (I.NE.O) PRINT = .FALSE. 
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BEGIN ITERATIONS 
SAVE X VALUES, THEN GET F VALUES 
NP= N + 14 
DO 100 IT= 1,MAXIT 
DO 10 IVBL= 1,N 
XSAVE(IVBL)= X(IVBL) 
10 CONT INUE 
CALL FCN(MCASE, EZ, LENGTH) 
TEST F VALUES AND SAVE THEM 
ITEST= O 
DO 20 IFCN= 1,N 
IF (ABS(F(IFCN)).GT.FTOL) ITEST= ITEST + 1 
FSAVE(IFCN)= F(IFCN) 
20 CONT INUE 
PRINT CURRENT VALUES IF PRINT IS .TRUE. 
IF(.NOT.PRINT) GO TO 30 
WRITE(6,1000)IT, X 
4000 FORMAT(/’ AFTER ITER. ’,13, 
+ ‘ X AND F VALUES ARE’ ,/,10E13.6) 
WRITE(6, 1001) F 
1001 FORMAT(/, 10E613.6) 
SEE IF FTOL IS MET. IF NOT, CONTINUE. IF SO, RETURN. 
30 IF(ITEST.NE.O) GO TO 35 
RETURN 
THIS DOUBLE LOOP COMPUTES THE PARTIAL DERIVATIVES OF EACH 
FUNCTION FOR EACH VARIABLE AND STORES THEM IN A COEFFICIENT 
ARRAY 
35 DO 50 JCOL=1,N 
X(JUCOL)= XSAVE(JCOL) + DELTA(UCOL) 
CALL FCN(MCASE, EZ, LENGTH) 
DO 40 IROW=1,N 
A(IROW,UCOL)= (F(IROW) - FSAVE(IROW))/DELTA(JCOL) 
40 CONT INUE 
RESET X VALUES FOR NEXT COLUMN OF PARTIALS 
X(JUCOL)= XSAVE(UJCOL) 
50 CONT INUE 
PUT NEGATIVE VALUES OF F AS R.H.S. AND CALL "ELIM" 


DO 60 IROW= 1,N 
A(IROW,NP)= -FSAVE(IROW) 
60 CONTINUE 
CALL ELIM(A, N, NP, 2) 
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Cc 
Cc 
Cc 


C 


CHECK IF COEFFICIENT MATRIX IS NOT TOO ILL-CONDITIONED 


DO 70 IROW= 1.N 


IF(ABS(A(TROW,IROW)).LE.1.£-5) GO To 998 


70 CONT INUE 


APPLY THE CORRECTIONS TO THE X VALUES, ALSO SEE IF XTOL 


IS MET. 


ITEST= O 
DO 80 IVBL=1,N 


X(IVBL)= XSAVE(IVBL) + A(IVBL,NP) 
IF(ABS(A(IVBL,NP)).GT.XTOL) ITEST= ITEST +1 


80 CONTINUE 


IF XTOL IS MET, PRINT LAST VALUES AND RETURN, ELSE DO 


ANOTHER ITERATION 


IF(ITEST.EQ.0O) GO 
100 CONTINUE 


TO 997 


MAXIT ITERATIONS HAVE BEEN DONE. 


RETURN 


XTOL IS MET. PRINT LAST VALUES. 


997 IF(.NOT.PRINT) GO 
WRITE (6, 1002)IT, 
1002 FORMAT(/,’ AFTER 


110 RETURN 


TO 110 
x 


IMERE © ,.135 
+ ’ X VALUES (MEETING XTOL) ARE’,/,10F13.5) 


PARTIALS FORM A NEARLY SINGULAR MATRIX. PRINT MESSAGE. 


998 WRITE(6, 1003) 


1003 FORMAT(/’ CANNOT SOLVE SYSTEM. MATRIX NEARLY SINGULAR 


STOP 


C NUMBER OF EQUATIONS IS INVALID. PRINT MESSAGE. 


Cc 


999 WRITE(6,1004)N 


1004 FORMAT(/’NUMBER OF EQNS PASSED TO NLSYST IS INVALID.’, 


+’ MUST BE 2<N<3. 
STOP 
END 


VALUE WAS 
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SUBROUTINE ELIM(AB, N, NP, NDIM) 


THIS SUBROUTINE SOLVES A SET OF LINEAR EQUATIONS. 

THE GAUSS ELIMINATION METHOD IS USED, WITH PARTIAL PIVOTING. 
MULTIPLE R.H.S. ARE PERMITTED, THEY SHOULD BE SUPPLIED 

AS COLUMNS THAT AUGMENT THE COEFFICIENT MATRIX. 

PARAMATERS ARE - 


AB COEFFICIENT MATRIX AUGMENTED WITH R.H.S. VECTORS 
N NUMBER OF EQUATIONS 
NP TOTAL NUMBER OF COLUMNS IN THE AUGMENTED MATRIX 


NDIM FIRST DIMENSION OF MATRIX AB IN THE CALLING PROGRAM. 
THE SOLUTION VECTOR(S) ARE RETURNED IN THE AUGMENTATION COLUMNS 
OF AB. 


DIMENSION AB(NDIM,NP) 
BEGIN THE REDUCTION 


NMi= N - 1 
0O 35 I= 1,NM1 


FIND THE ROW NUMBER OF THE PIVOT ROW. INTERCHANGE ROWS TO PUT 
THE PIVOT ELEMENT ON THE DIAGONAL 


IPVT= I 
IPi= I + 1 
DO 10 U= IP1, N 
IF(ABS(AB(IPVT,I)).LT.ABS(AB(JU,1I))) IPVT = JU 
40 CONTINUE 


CHECK IF THE PIVOT ELEMENT IS NOT TOO SMALL. IF SO PRINT 
A MESSAGE AND RETURN. 


IF ((ABS(AB(IPVT,I)).LT..00001)) GO TO 99 


INTERCHANGE, EXCEPT IF THE PIVOT ELEMENT IS ALREADY ON THE 
DIAGONAL . 


IF(IPVT.EQ.1I) GO TO 25 
DO 20 uvCOL= I,NP 
SAVE= AB(I,JCOL) 
AB(1I,JUCOL)= AB(IPVT,UJCOL) 
AB(IPVT,JCOL)= SAVE 
20 CONTINUE 
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REDUCE ALL ELEMENTS BELOW THE DIAGONAL IN THE I-TH ROW. 
CHECK FIRST TO SEE IF A ZERO ALREADY PRESENT. IF SO, 
CAN SKIP THE REDUCTION FOR THAT ROW 


25 DO 32 JROW= IP1,N 

IF(AB(JROW,I).EQ.0) GO TO 32 

RATIO= AB(JROW,1)/AB(I,1) 

DO 30 KCOL= IP1,NP 

AB(JUROW,KCOL)= AB(JROW,KCOL) - RATIO*AB(I,KCOL) 

30 CONT INUE 
32 CONT INUE 
35 CONTINUE 


CHECK A(N,N) FOR SIZE 
IF(ABS(AB(N,N)).LT..00001) GO TO 99 
BACK SUBSTITUTION 


NP1= N + 1 
DO 50 KCOL= NP1,NP 
AB(N,KCOL)= AB(N,KCOL)/AB(N,N) 
DO 45 J=2,N 
NVBL= NP1 - J 
L= NVBL + 1 
VALUE= AB(NVBL,KCOL) 
DO 40 K=L.N 
VALUE= VALUE - AB(NVBL,K)*AB(K,KCOL) 
40 CONTINUE 
AB(NVBL,KCOL)= VALUE/AB(NVBL ,NVBL) 
45 CONT INUE 
50 CONTINUE 
RETURN 


MESSAGE FOR A NEAR SINGULAR MATRIX 


99 WRITE(6, 100) 
100 FORMAT(/’SOLUTION NOT FEASIBLE. A NEAR ZERO PIVOT WAS 
+ ENCOUNTERED. ’ ) 
RETURN 
END 
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SUBROUTINE FCN (MCASE, EZ, LENGTH) 


THIS SUBROUTINE DEFINES ALL THE FUNCTIONS OF EQUILIBRIUM 
AND STRAIN COMPATIBILITY TO BE SATISFIED. 
PARAMETERS ARE - 
MCASE = 1 INDICATES THE ULTIMATE STAGE WHEN ANALYZING 
THE LOAD POINT SECTION 
= 2 INDICATES THE ULTIMATE STAGE WHEN ANALYZING 
THE CENTRE SUPPORT SECTION 
= 3 INDICATES THE UNCRACKED SECTION ANALYSIS 
EZ STEEL TENDON ECCENTRICITY 
LENGTH INDICATES THE INDEX IN THE VECTOR OF STEEL 
STRESS "FS" AT A DISTANCE FROM A CRACK 


DIMENSION X(2), F(2) 
COMMON X, F, B, H, EXL, EXS, PRESF, AREAP, CSTR, ELMOD, 


+FCU, E1, E2L, E2S, EPSO, SIGPU, AREAC, EPSU, COMPAT,DELTA(2), 


+ CTEN, ETEN, E2X, CRMOM, FS(60), STRAIN(60), RESULT, 
+ FRACL, SPAN, ISECT, BW, HF, R2, FSB, CURVCT, ASR, ETOP, 
+ RDFACT, RATMOD, ALPHA, BETA, GAMMA, AXF, EXTMOM, TCi, 
* ANGLEA, ANGLEB 


START THE SELECTION OF CASES 


ERCISECT EO. 0}° GOTO 20 
GO TO (11,22,33),MCASE 
20 GO TO (1,2,3),MCASE 


FOR MCASE = 1 AND 2, 
X(1) : LOCATION OF NEUTRAL AXIS 


X(2) : STEEL STRESS 
F(1) : EQUATION OF EQUILIBRIUM 
F(2) : STRAIN COMPATIBILITY EQUATION 
FOR MCASE = 3, 
X(1) : STRAIN OF EXTREME CONCRETE FIBRE IN COMPRESSION 
X(2) : STEEL STRESS 
F(1) : EQUATION OF EQUILIBRIUM 
F(2) : STRAIN COMPATIBILITY EQUATION 
AN : LOCATION OF NEUTRAL-AXIS 


RECTANGULAR SECTION 


14 F(1) = AREAP * X(2) * 1000. /(FCU*B) - X(1) 
F(2) = E1 + E2L + EPSU*(EXL+H/2.-X(1))/X(1) - PREL(X(2)) 
RETURN 

2 F(1) = AREAP * X(2) * 1000. /(FCU*B) - X(1) 
F(2) = E1 + E2S + EPSU*(EXS+H/2.-X(1))/X(1) - PREL(X(2)) 
RETURN 

3 AN = X(1)*H/(ETEN+X(1)) 
F(1) = ABS(COMP(AN,X(1))) - O.5*B*CTEN*(H-AN) + 

- AREAP*X(2)* 1000. 
F(2) = E1 + E2X + (ETEN+X(1))*(ABS(EZ)+H/2.)/H - PREL(X(2)) 
RETURN 
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I-SHAPED SECTION 


11 


22 


33 


+ 


ASR : AREA OF CONCRETE WITH THE WEB WIDTH TAKEN UP TO THE 


CROSS-SECTION TOP 


IF(X(2)*1000. *AREAP/(FCU*B) .LE. HF) GO TO 1 

ASR = AREAP - FCU*(B-BW)*HF/(X(2)*1000. ) 

F(1i) = 1.4*ASR*X(2)*1000./(BW*CSTR) - X(1) 

F(2) = E1 +E2L + EPSU*(EXL+H/2.-X(1))/X(1) - PREL(X(2)) 
RETURN 

IF(X(2)*1000.*AREAP/(FCU*B) .LE. HF) GO TO 2 

ASR = AREAP - FCU*(B-BW) *HF/(X(2)*1000. ) 

F(1) = 1.4*ASR*X(2)*1000./(BW*CSTR) - X(1) 

F(2) = E1 + E2S + EPSU*(EXS+H/2.-X(1))/X(1) - PREL(X(2)) 


AN = X(1)*H/(ETEN+X(1)) 

IF(AN .LE. HF) GO TO 3 

F(1) = ABS(COMP(AN,X(1))) - (BW*CTEN*(H-AN)/2. + 
(B-BW) *CTEN*HF) - AREAP*X(2)*1000. 


F(2) = E1 + E2X + (ETEN+X(1))*(ABS(EZ)+H/2.)/H - PREL(X(2)) 


RETURN 


END 
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SUBROUTINE AFCRAK(Y, EZ, KASE, LENGTH) 


SUBROUTINE FOR ROOT FINDING USING NEWTON’S METHOD. 


IT IS USED SPECIFICALLY FOR THE GENERAL CRACKED SECTION 


ANALYSIS. 
PARAMETERS ARE - 
Yi LOCATION OF THE NEUTRAL AXIS 
EZ STEEL TENDON ECCENTRICITY 
KASE = 1 INDICATES THE CONDITIONS JUST AFTER CRACKING 


AND AT A CRACK 


= 2 INDICATES THE CONDITIONS AT SOME DISTANCE FROM 


A CRACK 
LENGHT AS DEFINED IN SUBROUTINE "FCN" 


DIMENSION X(2), F(2) 


COMMON X, F, B, H, EXL, EXS, PRESF, AREAP, CSTR, 


ELMOD, 


+FCU bles EZ EQS ..ERSO.8 SIGPU, AREAC, EPSU, COMPAT,DELTA(2), 


+ CTEN, ETEN, E2X, CRMOM, FS(60), STRAIN(60), 


+++ 


ANGLEA, ANGLEB 


FRACL, SPAN, ISECT, BW, HF, R2, FSB, CURVCT, ASR, 
RDFACT, RATMOD, ALPHA, BETA, GAMMA, AXF, EXTMOM, TC1, 


RESULT, 
ETOP, 


START AN INITIAL VALUE AND CHECK IF EQUATION OF NEUTRAL AXIS 


IS SATISFIED 


IF(ISECT .EQ. O) BW=B 
FAXNT = AXNT(Y, EZ, KASE, LENGTH) 
DO..20, UC, =..1,,40 
DELY = FAXNT/DAXNT(Y, EZ, KASE, LENGTH) 
Y = Y - DELY 
FAXNT = AXNT(Y, EZ, KASE, LENGTH) 
IF(ABS(DELY) .LE. .0001) GO TO 70 
IF(ABS(FAXNT) .LE. .00001) GO TO 70 

20 CONTINUE 
WRITE (6,100) JC 


400 FORMAT (/’ TOL. NOT MET AFTER ’,15,’ ITERATIONS.’ ) 


STOP 
70 IF(Y .LT. O.) GO TO 30 
EVALUATION OF STRESSES 
FC : STRESS IN CONCRETE EXTREME FIBRE 
X(1) : CORRESPONDING STRAIN IN CONCRETE 
X(2) : STRESS IN PRESTRESSING STEEL 
FC = EXTMOM*Y/(GAMMA-BETA*Y-BW*Y**3/6. ) 
X(1) = FC/(57000. *SQRT(CSTR)) 
X(2) = RATMOD*(ABS(EZ)+H/2.-Y)*FC/Y + AXF/AREAP 
RETURN 
NEGATIVE NEUTRAL AXIS. PRINT MESSAGE. 
30 WRITE (6,204) Y 
204 FORMAT (/’ NEGATIVE N.A. = ’,€£12.5) 
STOP 
END 
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FUNCTION AXNT(Y, EZ, KASE, LENGTH) 


FUNCTION TO DEFINE THE EQUATION OF NEUTRAL AXIS 

IN THE GENERAL CRACKED SECTION ANALYSIS TO BE USED 
IN THE SUBROUTINE "AFCRAK" 

PARAMETERS ARE DEFINED IN "AFCRAK" 


DIMENSION X(2), F(2) 

COMMON X, F, B, H, EXL, EXS, PRESF, AREAP, CSTR, ELMOD, 

+FCU, E1, E2L, E2S, EPSO, SIGPU, AREAC, EPSU, COMPAT,DELTA(2), 

+ CTEN, ETEN, E2X, CRMOM, FS(60), STRAIN(60), RESULT, 

FRACL, SPAN, ISECT, BW, HF, R2, FSB, CURVCT, ASR, ETOP, 
RDFACT, RATMOD, ALPHA, BETA, GAMMA, AXF, EXTMOM, TCi, 
ANGLEA, ANGLEB 


+++ 


IF(ISECT .EQ. O) BW=B 
IF (KASE .EQ. 2) GO TO 2 


CONDITIONS JUST AT A CRACK 


AXNT = BW*AXF*Y**3/6. + BW*EXTMOM*Y**2/2. 
+ + (BETA*AXF + ALPHA*EXTMOM) *Y 
+ - (GAMMA*AXF + BETA*EXTMOM) 

RETURN 


CONDITIONS AT SOME DISTANCE FROM A CRACK 
SO THAT THE STEEL STRESS AND CONCRETE STRAIN CAN BE EVALUATED 
IN ORDER TO CALCULATE THE AVERAGE CURVATURE 


2 DEN = ABS(EZ)+H/2.-Y 
CHECK IF THE DENOMINATOR IS NEGATIVE. 


IF (DEN .EQ. 0.) DEN=DEN+.05 
TC1 = (BW*ty + (2*Y-HF) *HF*(B-BW)/Y)*FS(LENGTH) *Y 
+ /(2.*RATMOD*DEN) 
+ - AREAP*FS(LENGTH) - AXF/1000. 
AXNT = -AREAP*FS(LENGTH) *(ABS(EZ)+H/2.+ETOP) 
- TC1*(ETOP+(H+2.*Y)/3.) 
+ (.5*(Y/3.+ETOP) *BW*Y 
+ HF*(Y-HF)*(B-BW) *(HF/2.+ETOP)/Y 
+ .5*HF*HF*(B-BW)*(HF/3.+ETOP)/Y) 
* FS(LENGTH) *Y/(RATMOD*DEN) 


++ tet 


RETURN 
END 
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FUNCTION DAXNT(Y, EZ, KASE, LENGTH) 


FUNCTION DEFINING THE DERIVATIVE OF THE EQUATION OF 
NEUTRAL AXIS FOR USE IN THE SUBROUTINE "AFCRAK". 


THE 


DERIVATIVE OF THE EQUATION IS REQUIRED BY NEWTON’S 


ITERATIVE PROCEDURE (SLOPE). 


DIMENSION X(2), F(2) 
COMMON X, F, B, H, EXL, EXS, PRESF, AREAP, CSTR, ELMOD, 


+ 
+ 


+++ 


+t eet tet 


FCU, E1, E2L, E2S, EPSO, SIGPU, AREAC, EPSU, COMPAT,DELTA(2), 
CTEN, ETEN, E2X, CRMOM, FS(60), STRAIN(6O), RESULT, 
FRACL, SPAN, ISECT, BW, HF, R2, FSB, CURVCT, ASR, ETOP, 
RDFACT, RATMOD, ALPHA, BETA, GAMMA, AXF, EXTMOM, TC1, 
ANGLEA, ANGLEB 


IF (ISECT .EQ. O) BW=B 
IF (KASE .EQ. 2) GO TO 2 
DAXNT =BW*AXF*Y**2/2. +BW*EXTMOM*Y +(BETA*AXF +ALPHA*EXTMOM) 
RETURN 
DEN = ABS(EZ)+H/2.-Y 
IF (DEN .EQ. O.) DEN=DEN+.05 
TC = (BW*Y + (2.*Y-HF)*HF*(B-BW)/Y)*FS(LENGTH) *(ABS(EZ)+H/2. ) 
/(2.*RATMOD*DEN**2) 
+(BW +(B-BW) *HF*HF/Y**2)*FS(LENGTH) *Y/(2.*RATMOD*DEN) 
DAXNT = TC1*2./3. + (ETOP+(2.*Y+H)/3.)*TC 
- (.5*(Y/3.+ETOP) *BW*Y 
+ HF*(Y-HF)*(B-BW)*(HF/2.+ETOP)/Y 
+ .5*HF*HF*(B-BW)*(HF/3.+ETOP)/Y) 
*FS (LENGTH) *(ABS(EZ)+H/2.)/(RATMOD*DEN**2) 
- (.5*(BWtY/3.+ (Y/3.+ETOP) *BW) 
+HF *HF *(B-BW) *(HF/2.+ETOP)/Y**2 
- .S*HF *HF *(B-BW) *(HF/3.+ETOP)/Y**2) 
*FS(LENGTH) *Y/(RATMOD*DEN) 
RETURN 
END 
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FUNCTION PREL(SIG) 


FUNCTION DEFINING THE PRESTRESSING STEEL STRESS-STRAIN 
RELATIONSHIP FOR USE IN THE STRAIN COMPATIBILITY ANALYSIS 
AT THE ULTIMATE AND VARIOUS LOADING STAGES. 
THE STRESS-STRAIN CURVE HAS BEEN EXPRESSED IN AN 
EXPONENTIAL-TYPE EQUATION. 
PARAMETERS ARE - 

SIG STRESS OF STEEL 

THE VALUE OF THE CORRESPONDING STRAIN IS RETURNED. 


DIMENSION X(2), F(2) 

COMMON X, F, B, H, EXL, EXS, PRESF, AREAP, CSTR, ELMOD, 

+FCU, E1, E2L, E2S, EPSO, SIGPU, AREAC, EPSU, COMPAT,DELTA(2), 

CTEN, ETEN, E2X, CRMOM, FS(60), STRAIN(60), RESULT, 

FRACL, SPAN, ISECT, BW, HF, R2, FSB, CURVCT, ASR, ETOP, 
RDFACT, RATMOD, ALPHA, BETA, GAMMA, AXF, EXTMOM, TC1, 
ANGLEA, ANGLEB 


+ 


+ + + 


TwO TYPES ARE CONSIDERED: 
250 AND 270-KSI HIGH GRADE STEEL 


IF(SIGPU - 270.) 10,20,20 

10 DESIG pdgsadnt, 1,2 

2 PREL= SIG*1000./ELMOD + 2.5*((SIG-195.)**3)/10.**7 
RETURN 

20° IF(SIG -g210cdgd; 193 

3 PREL = SIG*1000./ELMOD + 2.*((SIG-210.)**3)/10. **7 
RETURN 


LINEAR RANGE 
4 PREL= SIG * 1000./ELMOD 


RETURN 
END 
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FUNCTION CMOM(C,EPS) 


FUNCTION CALCULATING THE CONTRIBUTION, TO THE TOTAL 
MOMENT, OF THE COMPRESSED CONCRETE AREA IN A SECTION. 
PARAMETERS ARE - 

Cc LOCATION OF THE NEUTRAL AXIS 

EPS STRAIN OF EXTREME CONCRETE FIBRE IN COMPRESSION 
THE FUNCTION TAKES INTO CONSIDERATION THE STRESS-STRAIN 
RELATIONSHIPS OF THE CONCRETE. 


DIMENSION X(2), F(2) 
COMMON X, F, B, H, EXL, EXS, PRESF, AREAP, CSTR, ELMOD, 
+FCU, E1, E2L, E2S, EPSO, SIGPU, AREAC, EPSU, COMPAT,DELTA(2), 
+ CTEN, ETEN, E2X, CRMOM, FS(60), STRAIN(60O), RESULT, 
FRACL, SPAN, ISECT, BW, HF, R2, FSB, CURVCT, ASR, ETOP, 
RDFACT, RATMOD, ALPHA, BETA, GAMMA, AXF, EXTMOM, TC1, 
ANGLEA, ANGLEB 


+++ 


IF (C .LT. O.) GO TO 50 

IF(EPS .GT. EPSU) GO TO 3 

FF = C*EPSO/EPS 

G = C-HF 

IF((ISECT .EQ. 0) .OR. (G .LE. O.)) GO TO 40 
IF(EPS .GT. EPSO) GO TO 20 


I-SHAPED SECTION. EPSO IS THE STRAIN AT WHICH THERE IS 
A CHANGE OF STRESS-STRAIN RELATIONSHIP. 
CASE WHEN THE APPLIED STRAIN IS BELOW EPSO 


CMOM = (B-BW)*CSTR*G**3*(8./3. - G/FF)/(4.*FF) + 
+ B*HF*CSTR*(HF*HF/3. + (C+G)**2 - 

+ ((C+G)**3 + (C+G)*HF**2)/(2.*FF))/(2. *FF) 
RETURN 


CASE WHEN THE APPLIED STRAIN IS HIGHER THAN EPSO 
201 (EF-G)27,21,23 
EPSO IS IN THE WEB 


21 CMOM = 5.*CSTR*(B-BW) *FF*FF/12.+ 
+ 0.5*CSTR*(B-BW) *(G-FF) *(G+FF-RDFACT* 
+(EPS*((G-FF) **2/3.+(GHFF) **2)/(2.*C)-EPSO*(G+FF))/(EPSU-EPSO) ) 
+ + 0.5*CSTR*B*HF * (C+G-RDFACT* 
+ (EPS*(HF**2/3.+(C+G)**2)/(2.*C)-EPSO*(C+G) )/(EPSU-EPSO) ) 


RETURN 
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EPSO IS IN THE FLANGE 


23 P = EPSO/EPS 
CMOM = (B-BW) *CSTR*G**3*(8./3.-G/FF)/(4.*FF) + 
. 5*B*CSTR*(FF-G) *(((FF+G) **2+(FF-G)**2/3.)/FF 
- ((FF+G) **3+(FF+G) *(FF-G) **2)/(4. *FF**2) ) 
+ .5*B*CSTR*C**2*(1.-P)*(1.+P- RDFACT* 
(EPS*((1.-P)**2/3.+(1.+P)**2)/2.-EPSO*(1.+P))/(EPSU-EPSO) ) 
RETURN 


++ ++ 


RECTANGULAR CROSS-SECTION 


40 IF(EPS .GT. EPSO) GO TO 2 
FF = EPS/EPSO 


CMOM l= Bet S*C)**21*CSTR*CFE*S./ 3.07 Iressa) 
RETURN 
2 FF = EPSO/EPS 
G = RDFACT/(EPSU-EPSO) 
CMOM = CSTR*S.*B*((C*FF)**2)/12. + B*((.5*C)**2) *CSTR*(1.-FF)* 
& COT SFE) (24 -Ge(ERS-EPSO)).-/ (C1 nehF)t*2)AERS*G/35) 
RETURN 


NEGATIVE VALUE OF NEUTRAL AXIS LOCATION. PRINT MESSAGE. 


5O WRITE (6,300) C 
300 FORMAT (/’ NEG. N.A. = ’,E£12.5,’ IN SECT. CMOM’ ) 
STOP 


THE ULTIMATE CONCRETE STRAIN HAS BEEN EXCEEDED. 
PRINT MESSAGE. 


3 WRITE(6,500) EPS 
500 FORMAT(//15X,F13.5,’ ,ULT CONC STRAIN EXCEEDED, ’ 
& /26X,’END OF COMPILING IN SECT CMOM.’) 
STOP 
END 
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FUNCTION COMP(C,EPS) 


FUNCTION CALCULATING THE CONTRIBUTION, TO THE TOTAL FORCE, 
OF THE CONCRETE AREA IN COMPRESSION. 
PARAMETERS ARE - 

Cc LOCATION OF THE NEUTRAL AXIS 

EPS STRAIN OF EXTREME CONCRETE FIBRE IN COMPRESSION 
THE FUNCTION TAKES INTO CONSIDERATION THE STRESS-STRAIN 
RELATIONSHIPS OF THE CONCRETE. 


DIMENSION X(2), F(2) 
COMMON X, F, B, H, EXL, EXS, PRESF, AREAP, CSTR, ELMOD, 


+FCU, E1, E2L, E2S, EPSO, SIGPU, AREAC, EPSU, COMPAT,DELTA(2), 


+ CTEN, ETEN, E2X, CRMOM, FS(60), STRAIN(60), RESULT, 
+ FRACL, SPAN, ISECT, BW, HF, R2, FSB, CURVCT, ASR, ETOP, 
+ RDFACT, RATMOD, ALPHA, BETA, GAMMA, AXF, EXTMOM, TCi1, 
+ ANGLEA, ANGLEB 


IF, (C. vl T= .09) GO*TOs50 
IF(EPS .GT.EPSU) GO TO 3 


I-SHAPED CROSS-SECTION. 


FF = C*EPSO/EPS 

G = C-HF 

IF((ISECT .EQ. 0) .OR. (G .LE. 0.)) GO TO 40 
IF(EPS .GT. EPSO) GO TO 20 


CASE WHEN APPLIED STRAIN IS BELOW EPSO 


COMP = (B-BW)*G*CSTR*(G/FF-(G/FF)**2/3. ) 
+ + Ries TRAHE*((C#G) /FF-((C#G) ##24HF#*2/3.)/(4.*FF**2)) 


RETURN 
20 IF(FF-G) 21, 21, 23 


EPSO IN THE WEB 


21 COMP= 2.*(B-BW)*FF*CSTR/3. + (B-BW) * 
+ (G-FF)*CSTR*(1.-RDFACT*(EPS*(G+FF)/C-EPSO) /(EPSU-EPSO) ) 
+o + B*HF*CSTR*(1.-RDFACT*(EPS*(C+G)/C-EPSO)/(EPSU-EPSO)) 


RETURN 


Tz 


Pet 


08s 


4 


0 Sas BIG ee aie es OTIS 


aie 


L = she aa 
20g? JATOT SMP (OY \MOy BR ow TA 


- 
vthe debian Sut 30. HOTT 

WOTZZAAINGD WI BRAal4 ar VOMFRTR? YO Mb 
WiaeTé-e2aave INT MOTT o ofut zaMat & 


y, e 
# 
< * 


: . ee sv Bor 

Oma RIBS. qn RA ae a. m2. a Ae - 
AT-236 HaFOD vueds oatehe 2 73. os * es .t ak 

aes toa braE (28% - eke _XS3 Wate ,MaT a a 
a0Ts aah “TOMRAGS 42) 62 aH (WR .TORST” Mage + 

2 MOMISS TER SMA eae Aaa GONTAR ,TORRGR + 
“wasova [3 iowa 4 | 


: o@ BY 68 (20 . ts. 3) 
an be “8 OF Gp CuBSS. TH. Raicaatii 


; Gb OT OF b4.0" Blo OY) ae. Be on 1388 x4 ee 
Pt ae C3 or uo Salt at, ps 


i + 


! Gae2 Wound ie winare Sh Sk a 
ae  stowtitheistenee + (yi = aNcD 
pes 3 ~ 


~ . Te Mao, ‘ 


5-8) + 
14 rtoedt he xt 7 
janes ays Sa ahah ih 


i 
oar a 
~ 440249 


aan 


aang0 


aqgagan qaaqgaqgn 


aqaag0 


V2 


EPSO IN THE FLANGE 


23 COMP = (B-BW)*G*CSTR*(G/FF-(G/FF)**2/3. ) 
+ + B*CSTR*(FF-G)*(1.+G/FF-(1.+G/FF+(G/FF)**2)/3.) 
+ + .5S*B*CSTR*C*(1.-EPSO/EPS) *(2.-RDFACT*(EPS-3.*EPSO)/ 
+ (2.*(EPSU - EPSO))) 
RETURN 


RECTANGULAR CROSS-SECTION 


40 IF(EPS.GT.EPSO) GO TO 2 

G= EPS/EPSO 
COMP= B*C*CSTR*(G-G*G/3. ) 
RETURN 

2 G= EPSO/EPS 
COMP= B*C*CSTR*(G*2./3. 

& + .5*(1.-G)*(2.-RDFACT*(EPS-EPSO)/(EPSU-EPSO) ) ) 
RETURN 


NEGATIVE VALUE OF LOCATION OF NEUTRAL AXIS. 
PRINT MESSAGE 


50 WRITE (6,300) C 
300 FORMAT (/’ NEG. N.A. = ’,£12.5,’ IN SECT. COMP’) 
STOP 


ULTIMATE CONCRETE STRAIN HAS BEEN EXCEEDED. 
PRINT MESSAGE. 


3 WRITE(6,501)EPS 
501 FORMAT(//15X,F13.5,’ ,ULT CONC STRAIN EXCEEDED, ’ 
& /26X,’END OF COMPILING IN SECT COMP. ’) 
STOP 
END 
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SUBROUTINE ROMB(IBOUND, EZ, DISTM) 


SUBROUTINE FOR ROMBERG INTEGRATION. PROGRAM BEGINS WITH 
TRAPEZOIDAL INTEGRATION WITH 10 SUBINTERVALS. INTERVALS 
ARE THEN HALVED AND RESULTS ARE EXTRAPOLATED UP TO 
FOURTH ORDER. 
MAXIMUM NUMBER OF SUBINTERVALS USED IN PROGRAM IS 160. 
PARAMETERS ARE - 

IBOUND NUMBER OF KNOWN VALUE INTEGRATION POINTS 

EZ PRESTRESSING STEEL ECCENTRICITY 

DISTM VALUE OF MOMENT 

RESULT RESULT OF INTEGRATION. RETURNED. 

TRAP DOUBLY SUBSCRIPTED ARRAY THAT HOLDS INTERMEDIATE 

VALUES FOR COMPARISONS AND EXTRAPOLATIONS 

OQ WHEN NON-CONVERGENT 
1 MEANS ALL OK. 


DIMENSION X(2), F(2) 
COMMON X, F, B, H, EXL, EXS, PRESF, AREAP, CSTR, ELMOD, 
+FCU, E1, E2L, E2S, EPSO, SIGPU, AREAC, EPSU, COMPAT,DELTA(2), 
+ CTEN, ETEN, E2X, CRMOM, FS(60), STRAIN(60), RESULT, 
FRACL, SPAN, ISECT, BW, HF, R2, FSB, CURVCT, ASR, ETOP, 
RDFACT, RATMOD, ALPHA, BETA, GAMMA, AXF, EXTMOM, TCi, 
ANGLEA, ANGLEB 


+++ 


DIMENSION TRAP(5,5) 
SET FLAG AT 1 INITIALLY 
KFLG = 1 


COMPUTE FIRST INTEGRAL WITH 10 SUBINTEGRALS AND USING TRAP 
RULE 


KINT = IBOUND/10 


SUM = STRAIN(1) + STRAIN(IBOUND) 
INT = O 
DO 10 I = 2,10 
INT = INT + KINT 

10 SUM = SUM + STRAIN(INT)*2. 
TRAP(1,1) = KINT/2*SUM 
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RECOMPUTE INTEGRAL WITH KINT HALVED, EXTRAPOLATE AND TEST. 
REPEAT UP TO 4 TIMES. 


DO 20 I = 1,4 
KINT = KINT/2 
INT = KINT 
K = 10*2**T 
DO 30 JU = 2,K,2 
SUM = SUM + STRAIN(INT)*2. 
INT = INT + KINT + KINT 
30 CONTINUE 
TRAP(1,1+1) = KINT/2*SUM 
DO 40 L = 1,I 
TRAP(L+1,I1+1) = TRAP(L,I+1) + 1./(4.**L - 1.)* 
& (TRAP(L,I+1) - TRAP(L,I)) 
40 CONTINUE 
IF (ABS(TRAP(I+1,1+1) - TRAP(I,I+1)) - .01) 50,50,20 
20 CONTINUE 


IF TOLERANCE NOT MET AFTER 4 EXTRAPOLATIONS, PRINT 
MESSAGE. SET KFLG = O 

KFLG = O 

WRITE (6,200) 


200 FORMAT (/’TOLERANCE NOT MET. CALCULATED VALUES WERE ’) 


SOhI-= I +1 
IF (KFLG .EQ. O) STOP 


CALCULATION OF FINAL RESULT 


RESULT = (DISTM/ABS(DISTM))*TRAP (1,1) / 

+ (ABS(EZ)+H/2.) 
RETURN 
END 
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